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Abstract:  Because  of  the  observed  decrease  in  the  ice  cover  in  the  Beau¬ 
fort  and  Chukchi  Seas,  it  is  possible  that  spray  icing  of  rigs  used  for  oil  ex¬ 
ploration  and  drilling  may  be  more  frequent  and  possibly  more  severe  in 
the  coming  years  than  it  has  been  in  the  past.  In  this  report  we  describe  a 
model  for  sea  spray  icing  on  fixed  offshore  structures.  The  accretion  of 
small  sea  spray  droplets  onto  two-dimensional  structural  sections  and 
components  depends  on  the  liquid  water  content  of  the  spray  cloud,  as 
well  as  wind  speed,  droplet  diameter,  and  the  diameter  of  the  object.  The 
spray  cloud’s  liquid  water  content  is  obtained  from  the  flux  of  film,  jet,  and 
spume  droplets  from  the  ocean  surface  and  the  vertical  velocity  of  the 
droplets.  The  spray  droplet  flux  increases  dramatically  with  increasing 
wind  speed,  as  whitecaps  cover  more  of  the  ocean  surface  and  the  wind 
shears  droplets  off  the  wave  crests.  The  more  massive  larger  droplets  tend 
to  fall  out  of  the  spray  cloud;  as  a  result,  the  liquid  water  content  decreases 
with  height.  We  present  modeled  icing  rates  for  the  semi-submersible 
drilling  rigs  Ocean  Bounty  in  the  winter  of  1979-1980  and  the  Sedco  708 
in  January  1983.  These  results  are  compared  to  available  information  on 
the  icing  rate  or  the  ice  accumulation. 
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1  Introduction 

Because  of  the  observed  decrease  in  the  ice  cover  in  the  Beaufort  and 
Chukchi  Seas,  it  is  possible  that  spray  icing  of  rigs  used  for  oil  exploration 
and  drilling  may  be  more  frequent  and  possibly  more  severe  in  the  coming 
years  than  it  has  been  in  the  past.  While  a  complete  ice  cover  can  result  in 
substantial  ice  forces  on  offshore  structures,  waves  are  not  a  problem.  In 
waters  with  a  partial  ice  cover,  wave  heights  tend  to  be  small,  even  in  high 
winds.  However,  the  overall  decrease  in  extent  of  the  sea  ice  cover  has  left 
open  water  in  parts  of  the  Beaufort  and  Chukchi  Seas  that  used  to  be  ice 
covered.  A  comparison  of  the  sea  ice  extent  in  September  1985  and  Sep¬ 
tember  2005  is  shown  in  Figure  1 

f  http://nsidc.org/data/seaice  index/archives/index.htmO.  Oil  exploration  platforms  and 
drill  rigs  operating  in  these  areas  now  experience  more  open  water  and 
higher  waves  and,  thus,  could  increasingly  suffer  from  sea  spray  icing.  It 
would  therefore  be  prudent  to  provide  an  icing  forecast  for  these  offshore 
structures  operating  in  open  water  during  subfreezing  temperatures,  simi¬ 
lar  to  that  now  provided  for  vessels  by  the  National  Oceanic  and  Atmos¬ 
pheric  Administration  (NOAA).  This  forecast  would  allow  the  rig  operators 
to  change  or  curtail  operations  and  to  prepare  to  operate  any  available  de¬ 
icing  equipment  on  crucial  portions  of  the  structure  and  gear. 


Figure  1.  Change  in  sea  ice  extent  from  September  1985  to  September  2005.  The  pink  line 
shows  the  median  ice  edge  in  September. 
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Platform  icing  versus  ship  icing 

Sea  spray  icing  of  fishing  and  crabbing  boats,  coastal  freighters,  Navy  and 
Coast  Guard  ships,  and  tenders  for  offshore  structures  in  Arctic  and  north¬ 
ern  waters  in  the  winter  is  a  well-known  hazard.  The  weight  of  ice  accreted 
in  the  superstructure  of  relatively  small  vessels  can  lead  them  to  lose  sta¬ 
bility.  In  these  conditions,  the  vessel  often  sinks  rapidly  with  the  loss  of  all 
hands.  On  larger  vessels,  in  the  same  conditions,  the  severity  of  super¬ 
structure  icing  is  less,  with  the  threshold  wind  speed  required  to  cause 
spray  icing  increasing  with  the  length  of  the  ship  (Table  5,  Overland  1990). 

NOAA  has  produced  an  ice  accretion  forecast  chart  on  a  daily  basis  since 
the  winter  of  1986-1987  using  the  Overland  et  al.  (1986)  algorithm  with  a 
predictor  PR  that  depends  on  wind  speed  V,  sea  temperature  Tw ,  and  air 

temperature  Ta : 


PR 


V(Tf-T') 

1  +  0A(TW-Tf)' 


(1) 


where  77  is  the  salinity-dependent  freezing  point  of  seawater  (-i.8°C  in  the 

north  Pacific).  The  predictor  is  related  to  icing  rate  (cm  Ir1)  on  20-  to  75- 
m-long  trawlers  that  are  headed  into  or  abeam  of  the  wind  and  not  at¬ 
tempting  to  ameliorate  the  accretion  of  ice  either  by  avoiding  open  seas  or 
by  reducing  speed.  The  numerator  in  Equation  1  represents  cooling  of  the 
spray  water  by  the  air,  and  the  terms  in  the  denominator  represent  the 
heat  of  fusion  of  the  spray  that  freezes  and  the  cooling  of  all  the  spray  wa¬ 
ter  to  the  freezing  point.  For  the  ice  accretion  forecasts,  air  temperature 
and  wind  speed  are  the  2-m  and  10-m  values,  respectively,  from  the  Global 
Forecast  System  model.  The  sea  temperature  is  obtained  from  a  satellite, 
ship,  and  buoy  sea-surface  temperature  analysis.  Forecast  ice  accretions 
are  mapped  as  light  (<0.7  cm  Ir1),  moderate  (0.7-2. o  cm  h1),  or  heavy  (>2 
cm  Ir1)  on  a  10  by  10  grid.  A  recent  example  is  shown  in  the  National  Cen¬ 
ter  for  Environmental  Prediction  12-hour  forecast  for  0000  on  November 
6,  2006,  in  Figure  2. 

Sea  spray  icing  on  stationary  offshore  structures,  including  oil  exploration 
platforms  and  drill  rigs,  is  significantly  different  from  sea  spray  icing  on 
ships.  Spray  is  generated  on  ships  by  heaving  and  pitching  as  the  ship  in¬ 
teracts  with  the  waves  it  is  moving  through.  How  individual  vessels  per¬ 
form  in  waves  depends  on  their  length,  freeboard,  stability,  and  hull 
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shape;  so  the  icing  severity  depends  on  those  characteristics  of  the  vessel 
as  well  as  on  the  environmental  conditions.  The  majority  of  the  spray  that 
causes  ship  icing  comes  from  the  impact  of  the  water  against  the  hull  as 
the  vessel  slams  into  the  waves.  On  the  other  hand,  offshore  structures  for 
exploration  and  drilling,  including  submersibles,  semi-submersibles,  and 
jack-up  drilling  vessels,  are  essentially  stationary,  either  sitting  directly  on 
the  bottom  or  anchored.  They  also  tend  to  be  relatively  open  at  the  water 
level  in  contrast  to  the  near-vertical  sides  of  ships.  Thus  the  spray  that  af¬ 
fects  these  offshore  structures  tends  to  be  generated  from  the  waves  them¬ 
selves  rather  than  from  the  structure-wave  interaction. 


Figure  2.  Ice  accretion  forecast  for  ships.  Units  are  cm  h1. 


The  goal  of  this  study  is  to  provide  algorithms  for  processing  weather  data 
to  determine  sea  spray  icing  severity  on  offshore  structures.  Input  to  the 
algorithms  is  data  from  offshore  platforms,  buoys,  or  stations  on  the  coast 
with  hourly  measurements  of  wind  speed  and  air  temperature.  For  coastal 
stations,  only  data  from  periods  with  onshore  winds  should  be  used  so  that 
over-ocean  conditions  are  represented.  Ultimately,  these  algorithms  could 
be  used  with  forecast  wind  speeds  and  temperatures  to  forecast  icing  se¬ 
verity  on  drill  rigs,  similar  to  the  vessel  icing  forecasts  that  are  currently 
available  (Figure  2).  With  such  a  forecast,  the  operators  of  drill  rig  and  oil 
exploration  platforms  could  prepare  to  deploy  anti-icing  and  deicing 
measures  and  curtail  operations  to  an  extent  determined  by  the  severity  of 
the  expected  icing  conditions.  Icing  is  a  safety  issue  for  workers  on  the  ves¬ 
sels  and  can  also  immobilize  exposed  mechanical  devices  and  thereby  limit 
both  communications  and  operations. 
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Organization 

In  Section  2  we  review  the  formulation  for  the  accretion  of  small  droplets 
on  two-dimensional  structural  sections  (e.g.,  cables,  guys,  angles,  chan¬ 
nels,  W-sections)  and  components  (e.g.,  handrails,  ladders,  antennas, 
winches).  Particular  attention  is  paid  to  the  collision  efficiency  of  the  spray 
droplets  and  the  effects  of  wind  speed,  droplet  diameter,  and  the  diameter 
of  the  component.  The  use  of  the  median  volume  diameter  to  characterize 
the  spray  droplet  cloud  is  discussed.  In  Section  3  we  introduce  the  formu¬ 
lation  for  the  flux  of  spray  droplets  from  the  ocean  surface.  The  physical 
mechanisms  of  spray  droplet  formation  are  briefly  described  and  refer¬ 
ences  with  more  detailed  discussions  are  provided.  The  empirical  relation¬ 
ships  between  significant  wave  height  and  wind  speed  as  well  as  whitecap 
generation  and  wind  speed  are  presented  and  their  limitations  are  dis¬ 
cussed.  The  calculation  of  the  spray  droplet  concentration  distribution  at 
the  sea  surface  from  the  spray  flux  is  discussed  and  the  expected  vertical 
variation  in  the  spray  droplet  profile  is  derived.  The  calculated  concentra¬ 
tion  distributions  at  half  the  significant  wave  height  are  compared  to 
measured  values  from  the  literature.  In  Section  4  we  present  model  results 
for  sea  spray  icing  on  the  Ocean  Bounty  semi-submersible  exploratory 
drilling  rig  in  the  winter  of  1979-1980  and  the  Sedco  708  semi- 
submersible  drilling  rig  in  January  1983.  These  results  are  compared  to 
available  information  on  the  icing  rate  or  the  ice  accumulation.  Limita¬ 
tions  of  the  current  state  of  knowledge  and  suggestions  for  further  studies 
are  discussed  in  Section  5. 
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2  Ice  Accretion 

Sea  spray  icing  on  fixed  offshore  platforms  occurs  when  wind-blown  water 
droplets  generated  from  whitecaps  on  the  ocean  surface  strike  a  portion  of 
the  platform  structure  and  freeze  (Figure  3).  We  expect  icing  to  be  more 
severe  at  lower  heights  above  the  ocean  surface,  both  because  the  droplets 
tend  to  evaporate  as  they  are  transported  from  their  source  and  because 
larger  droplets  tend  to  fall  out  of  the  spray  cloud  because  of  gravity.  On  the 
other  hand,  no  ice  is  expected  to  accrete  near  the  waterline,  where  the 
structure  is  warmed  by  the  water  as  waves  and  swell  wash  over  it.  Because 
the  spray  droplets  are  carried  by  the  wind,  more  ice  is  expected  to  accrete 
on  the  windward  side  of  the  platform  and  on  the  windward  side  of  struc¬ 
tures  and  components  on  the  platform.  Ice  will  also  accumulate  on  top  of 
horizontal  surfaces  as  water  flows  off  vertical  surfaces  or  drips  off  cables, 
handrails,  and  other  elevated  components  when  it  is  not  cold  and  windy 
enough  to  freeze  the  impacting  droplets  immediately. 


Figure  3.  Sea  spray  icing  on  the  Ocean  Bounty  in  the  winter  of  1979-1980  in  Lower  Cook  Inlet 

(photo  provided  by  MMS). 

The  spray  cloud  is  characterized  by  its  liquid  water  content  distribution 
(units  are  g  nr  3  j^nr1): 
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W 


4 

3 


^3P  wc{z,r) 


(l) 


where  pw  is  the  density  of  sea  water  (kg  nr3),  z  is  height  above  the  ocean 
surface  (m),  and  r  is  the  droplet  radius  (pm).  The  concentration  distribu¬ 
tion  of  the  spray  droplets  c(z,r)  is  the  number  of  droplets  per  cubic  meter 
of  air  per  micron  increment  in  radius  (nr 3  pm1)-  The  important  meteoro¬ 
logical  parameters  are  wind  speed  and  air  temperature.  The  most  severe 
icing  will  occur  in  events  with  high  winds  and  subfreezing  temperatures 
for  extended  periods. 

The  icing  rate  on  a  cylinder  is 


dTjz) 

dt 


U{z)D 


Pz 


'max 

I 


E(U ,r,D)w[z,r)dr , 


(2) 


where: 

T(z)  =  ice  thickness  (mm) 
t  =  time  (s) 

C7(z)  =  wind  speed  (m  S'1) 

D  =  cylinder  diameter  (m) 
p;  =  ice  density  (kg  nr3) 

E  {U,r,  D)  =  collision  efficiency  of  the  droplets  with  the  cylinder. 

The  collision  efficiency  is  a  number  between  o  and  l  that  is  the  portion  of 
the  mass  of  droplets  that  are  swept  out  by  the  cylinder  that  actually  hit  the 
cylinder.  It  results  from  a  balance  between  inertia,  which  tends  to  keep  the 
droplets  traveling  in  a  straight  line,  and  drag,  which  tends  to  make  the 
droplets  follow  the  wind  streamlines  around  the  cylinder.  Collision  effi¬ 
ciencies  for  cylinders  were  computed  on  an  analog  computer  and  provided 
in  tabular  graphical  form  in  Langmuir  and  Blodgett  (1946)  and  recom¬ 
puted  with  better  drag  formulations  on  a  digital  computer  and  provided  as 
convenient  formulas  by  Finstad  et  al.  (1988a). 

A  simplification  of  Equation  2  is  obtained  by  using  the  median  volume  ra¬ 
dius  ?’m,  which  is  the  radius  for  which  half  the  liquid  water  content  W(z)  is 
in  larger  droplets  and  half  is  in  smaller  droplets: 
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W  z 


'  M 

J  w[z,r)dr . 


(3) 


It  has  been  shown  for  cloud  droplet  distributions  (Finstad  et  al.  1988b) 
that 


r 

'  max 

f  E(U,r,D)w(z,r)dr  «  E(U,rM,D)W(z) ,  (4) 

rmin 

so  that  Equation  2  simplifies  to 

dT(z)  _U(z)DE(U,rM,D)W(z) 
dt  Pf 


The  range  of  drop  radii  from  sea  spray  (nominally  0.1-500  pm)  is  much 
larger  than  the  typical  range  in  clouds  (1-50  pm),  so  in  practice  Equations 
4  and  5  are  applied  over  segments  of  the  spray  droplet  distribution.  The 
liquid  water  content  distribution,  characterized  by  the  distribution  at  n 
specific  radii  r,-,  is  integrated  numerically  to  obtain  the  liquid  water  con¬ 
tent  using  the  trapezoid  rule.  The  median  volume  radius  rMj  for  that  por¬ 
tion  of  the  distribution  is 


'Mj 


=  rM  + 


WJ->  jl 

2  ,  1 

WJ~WM  V 

2 

0-0-1 


(6) 


where  we  use  the  plus  sign  when  Wj  >  Wj~i  and  the  minus  sign  when  Wj  < 
Wj-i.  The  collision  efficiency  is  then  determined  for  each  rMj,  and  the  icing 
rate  in  Equation  5  is  a  summation  over  the  n  droplet  radii: 


dTjz) 

dt 


U(z)D^E(U,r„J,D)WJ(z) 

3= 1 


P; 


(7) 


The  icing  rate  on  representative  structural  components  is  straightforward 
to  compute  given  the  concentration  distribution  of  droplets  in  the  spray 
cloud.  Because  the  droplets  cool  rapidly  from  the  sea  temperature  to  a 
temperature  lower  than  the  air  temperature,  a  heat  balance  calculation, 
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such  as  was  done  for  freezing  rain  in  Jones  (1996),  is  not  necessary.  This 
dry  accretion  (i.e.,  flux  limited)  assumption  (Jones  1998)  may  be  slightly 
conservative  at  air  temperatures  near  freezing,  causing  us  to  overestimate 
the  amount  of  ice  that  accretes.  However,  its  simplicity  more  than  com¬ 
pensates  for  these  occasional  overestimates. 

At  typical  ocean  salinities,  the  freezing  point  is  -i.8°C.  The  droplets  will 
tend  to  reject  brine  as  they  freeze,  leaving  a  salty  liquid  surface  layer  on 
the  frozen  droplet.  Except  in  relatively  warm  conditions  with  high  liquid 
water  contents,  the  droplets  freeze  individually  incorporating  small  air  in¬ 
clusions.  Thus  the  ice  accumulates  as  an  ice-water-air  matrix  with  a  den¬ 
sity  that  depends  on  the  weather  conditions  and  the  location  on  the  plat¬ 
form.  Ryerson  and  Gow  (2000)  measured  the  density  of  ice  accreted  on 
the  superstructure  of  the  Coast  Guard  cutter  Midgett  in  the  Bering  Sea  in 
1990.  The  density  of  17  samples  taken  from  various  locations  ranged  from 
690  to  920  kg  nr3,  with  a  median  of  850  kg  m  3.  They  also  summarize  den¬ 
sity  data  from  other  authors,  reporting  an  overall  range  from  620  to  940 
kg  m-3.  For  this  study  we  do  not  attempt  to  model  the  ice  density  but  as¬ 
sume  a  relatively  high  density  of  900  kg  m  3.  At  cold  temperatures  and 
relatively  low  spray  concentrations,  densities  are  likely  to  be  lower  than 
this.  In  these  conditions,  which  can  occur  with  moderate  whitecap  cover¬ 
age  at  low  heights  on  the  drill  rig  or  with  significant  whitecap  coverage 
higher  on  the  rig,  the  thickness  of  lower  density  ice  will  be  greater  than  the 
modeled  thickness.  However,  the  modeled  thickness  with  the  assumed 
density  correctly  represents  the  accreted  mass  of  ice. 

The  icing  rate  profile  in  Equation  7  is  calculated  throughout  the  spray  icing 
event  using  the  time-varying  measured  wind  speed  and  air  temperature 
and  modeled  droplet  concentration  density.  To  estimate  the  wind  speed 
profile,  we  use  the  standard  formulation  for  the  atmospheric  surface  layer, 
which  is  based  on  the  air-sea  temperature  difference  and  the  measured 
wind  speed  at  one  height  (e.g.,  Andreas  1998).  Andreas  et  al.  (2008)  de¬ 
scribe  the  specific  algorithm  that  we  use.  The  spray  generation  function, 
from  which  the  droplet  concentration  is  obtained,  increases  as  the  third 
power  of  wind  speed.  The  collision  efficiency  E  also  increases  with  wind 
speed.  Thus  the  rate  of  ice  accretion  in  Equation  7  increases  with  at  least 
the  fourth  power  of  wind  speed,  and  the  icing  rate  for  a  wind  speed  of  20 
m  s1  is  more  than  16  times  the  icing  rate  for  a  wind  speed  of  10  m  sA 
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3  Spray 

Sea  spray  fluxes 

An  icing  forecast  for  fixed  structures  is  based  on  the  potential  for  sea  spray 
generation  by  the  wind.  Sea  spray  droplets  are  produced  by  three  mecha¬ 
nisms  (Andreas  2002).  Wind  creates  breaking  waves;  these  waves  trap  air 
that  manifests  as  whitecaps  on  the  surface.  These  whitecaps  are  clouds  of 
bubbles;  when  the  bubbles  rise  to  the  surface  and  burst,  they  fling  into  the 
air  small  “film”  droplets  and  larger  “jet”  droplets  (Monahan  1986; 
Monahan  et  al.  1986).  The  small  film  droplets,  ranging  from  0.5  to  50  pm 
in  radius  are  produced  as  air  bubbles  in  a  whitecap  burst  and  the  bubble 
film  shatters  into  water  droplets.  Jet  droplets,  ranging  from  1  to  100  pm  in 
radius,  are  produced  as  the  collapsing  whitecap  bubble  shoots  up  a  jet  of 
water  that  breaks  into  droplets.  Spume  droplets,  also  called  spindrift,  are 
formed  by  the  wind’s  shearing  droplets  off  the  wave  crests;  these  range  in 
radius  from  20  to  500  pm.  The  peak  in  the  volume  flux  of  droplets— that 
is,  the  range  that  carries  most  of  the  airborne  water— occurs  in  the  spume 
size  range.  At  about  10  m  s_1,  the  winds  are  strong  enough  to  create  white- 
caps  covering  1%  of  the  sea  surface  and  tear  droplets  off  the  wave  crests. 
With  winds  of  20  m  s_1,  whitecaps  cover  10%  of  the  sea  surface.  Because 
spray  generation  follows  whitecap  coverage,  it  increases  at  roughly  the 
third  power  of  wind  speed  (Monahan  et  al.  1986). 

The  key  piece  of  information  for  this  study  is  the  rate  at  which  spray  drop¬ 
lets  are  produced.  This  spray  generation  function  is  denoted  dF/ dr,  and 
has  units  of  number  of  droplets  produced  per  square  meter  of  sea  surface 
area,  per  second,  per  micron  increment  in  droplet  radius  (nr2  s1  pm1)-  In 
the  field  of  marine  aerosol  science  there  has  been  quite  a  bit  of  controversy 
as  to  the  magnitude  and  form  of  the  spray  generation  function.  Significant 
advances  in  understanding  and  quantifying  sea  spray  generation  occurred 
as  a  result  of  the  HEXOS  (Humidity  Exchange  Over  the  Sea)  program  that 
began  in  the  early  1980s.  Andreas  (2002)  reviewed  about  a  dozen  pub¬ 
lished  sea  wave  spray  generation  functions.  He  used  known  relationships 
among  the  different  droplet  sizing  conventions  to  convert  each  of  the  pub¬ 
lished  spray  generation  functions  into  terms  of  the  spray  droplet  radius  at 
creation  at  the  level  of  the  sea  surface  r0.  Some  of  the  functions  had  to  be 
adjusted  from  measurement  heights  of  tens  to  hundreds  of  meters  above 
the  sea  surface.  He  made  these  adjustments  using  a  formulation  provided 
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by  Fairall  and  Larsen  (1984)  for  the  vertical  transport  of  droplets.  Andreas 
discarded  most  of  the  published  functions  because  they  violated  theoreti¬ 
cal  constraints  on  the  magnitude  of  the  function  and  its  wind  speed  de¬ 
pendence  and  concluded  that  the  spray  generation  function  that  Fairall  et 
al.  (1994)  developed  had  the  best  properties. 

The  Fairall  et  al.  function  estimates  the  rate  of  formation  of  spray  droplets 
with  r0  between  1.6  and  500  pm.  For  this  study  we  extended  it  to  even 
smaller  droplets  using  the  “bubbles-only”  function  reported  by  Monahan 
et  al.  (1986).  It  treats  film  droplets  as  small  as  0.5  pm  and  has  proven  use¬ 
ful  in  many  other  studies  of  marine  aerosols  (e.g.,  Burk  1984;  Stramska 
1987;  Gong  et  al.  1997).  Moreover,  because  it  has  the  same  wind  speed  de¬ 
pendence  as  the  Fairall  et  al.  (1994)  function,  it  was  easy  to  create  a  spray 
generation  function  that  treats  droplets  with  r0  from  0.5  to  500  pm  by 
simply  matching  the  Monahan  et  al.  and  Fairall  et  al.  functions  at  r0  =  2.0 
pm. 

Because  we  are  ultimately  interested  in  the  mass  flux  of  the  spray,  we  usu¬ 
ally  present  the  spray  generation  function  as  a  volume  flux 
(4  7r  r03  /  3)  dF  /  dr0 .  The  merged  spray  generation  volume  flux  is  plotted  in 

Figure  4  for  droplet  radii  r0  between  0.5  and  500  pm  for  four  wind  speeds. 
The  equilibrium  radius  for  droplets  in  air  of  80%  relative  humidity  rg 0  is 
roughly  half  of  r0;  rso  «  o.5r0.  The  figure  shows  that  the  rate  at  which  spray 
is  generated  increases  dramatically  with  wind  speed.  The  figure  also  sug¬ 
gests  that  droplets  with  radii  of  100  to  200  pm  carry  most  of  the  mass 
away  from  the  sea  surface.  These  droplets  are  almost  all  spume  droplets— 
those  ripped  right  off  the  wave  crests  by  the  wind. 

The  concentration  distribution  of  the  droplets  c(z,r)  in  the  air  depends  on 
the  height  above  the  water  surface  z  and  the  droplet  radius  r  and  can  be 
determined  from  the  concentration  at  the  ocean  surface.  Spray  droplets 
start  with  initial  radius  r0  and  with  the  same  temperature  as  the  surface 
ocean  water  Tw  but  begin  evolving  in  both  temperature  and  size  once  they 
are  ejected  into  the  air.  Figure  5  shows  the  evolution  of  a  droplet  with  100- 
pm  initial  radius  created  in  water  that  is  at  the  freezing  point  -i.8°C  and 
then  ejected  into  air  with  a  much  lower  temperature,  say,  -io°C.  These  re¬ 
sults  are  from  the  microphysical  spray  model  that  Andreas  (1989)  devel¬ 
oped  from  ideas  and  equations  in  Pruppacher  and  Klett  (1978). 
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Figure  5  makes  several  essential  points.  First,  this  100-pm  droplet  cools 
from  the  water  temperature  to  a  temperature  below  the  ambient  air  tem¬ 
perature  in  about  2  s.  The  smaller  droplets  cool  to  temperatures  below  the 
ambient  air  temperature  in  even  less  time.  In  typical  Arctic  conditions,  the 
spray  droplets  are  thus  supercooled  and  will  freeze  quickly  if  they  hit  a 
structure  with  a  surface  temperature  that  is  below  freezing. 

The  second  point  that  Figure  5  makes  is  that,  though  the  droplets  cool  very 
rapidly,  they  evaporate  much  more  slowly— the  evaporation  rate  is  about 
three  orders  of  magnitude  slower  than  the  cooling  rate  regardless  of  drop¬ 
let  size.  Furthermore,  because  the  droplets  are  saline,  they  do  not  evapo¬ 
rate  entirely  but,  rather,  reach  an  equilibrium  radius  that  is  no  smaller 
than  half  the  original  radius  for  typical  humidities  and  Arctic  Ocean  salini¬ 
ties.  For  example,  the  100-pm  droplet  in  Figure  5  retains  nearly  its  original 
radius  for  almost  100  s,  and  takes  2000  s  to  evaporate  to  a  60-pm-radius 
droplet.  For  this  study  we  assume  that  the  droplets  not  only  cool  to  the  air 
temperature  but  also  evaporate  to  their  equilibrium  size  at  80%  humidity 
before  being  blown  onto  the  superstructure  of  an  offshore  platform. 


Figure  4.  The  spray  generation  function  for  this  study  is  expressed  in  terms  of  a  volume  flux 
(4-rcr8o  / 3 )dF /  dr8o ,  where  rso  is  the  droplet  radius  in  equilibrium  at  a  reference  relative 

humidity  of  80%.  This  function  is  a  synthesis  of  functions  reported  by  Monahan  et  al.  (1986) 
and  Fairall  et  al.  (1994).  U10  is  the  wind  speed  at  a  height  of  10  m. 
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Time  Since  Formation  (s) 


Figure  5.  The  evolution  in  temperature  and  radius  of  a  droplet  with  100-pm  initial  radius  that 
formed  from  water  with  temperature  -1.8 °C  and  is  ejected  into  air  with  temperature  -10 °C 
and  relative  humidity  90%.  The  barometric  pressure  is  assumed  to  be  1000  mb,  and  the 

surface  salinity  is  34  psu. 

Significant  wave  height 

The  significant  wave  height  H{/?J  is  the  average  height  of  the  highest  one- 
third  of  the  waves.  The  height  of  wind  waves  increases  with  the  square  of 
the  wind  speed  and  also  depends  on  the  water  depth  Dw.  Andreas  and 
Wang  (2007)  determined  a  relationship  between  U10  (m  s1)  and  iTi/3  (m) 
based  on  buoy  data  in  the  North  Atlantic  off  the  east  coast  of  the  United 
States: 

Hi/3=c(dw)  for  ^10  <  4  m  s-1 

H!/3  =  a(Dw)uJ  +b(Dw)  for  U10  >  4  rn  s-1 

with 

,  ,  \ln(Dw/6) 

c(D  )  =  1.36 tanh  — - - -  (8; 

v  '  1.9 

a(Dw)  =  0.0134 tanh 

t-3 

b(Dw)  =  c(Dw)-i6a(Dw). 


The  measured  significant  wave  heights  are  scattered  about  the  values  cal¬ 
culated  from  Equation  8  because  of  variations  in  fetch,  swell,  the  stability 
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of  the  atmosphere,  and  the  duration  of  the  wind.  It  is  not  known  how  well 
this  parameterization  represents  conditions  off  the  Alaskan  coast. 

Whitecap  coverage 

The  availability  of  spray  droplets  is  determined  by  the  whitecap  cover  on 
the  surface  of  the  ocean.  For  wind  speeds  below  about  3  m  S’1  there  are  no 
whitecaps.  The  fraction  of  the  ocean  surface  that  is  covered  by  whitecaps 
fwc  is  estimated  from  the  wind  speed  at  10  m  using  a  relation  from 
Monahan  and  Muircheartaigh  (1980) 

fwc=3.84xw6U10™\  (9) 

Note  that  100%  whitecap  cover  is  obtained  for  Uw=39  m  s1.  The  data  un¬ 
derlying  Equation  9  is  for  relatively  warm  water  (20°C  <  Tw<  30°C),  and 
U10  <  20  m  S'1.  Those  authors  suggest  that  the  water  temperature  affects 
the  mean  lifetimes  of  whitecaps,  and  the  air-water  temperature  difference 
also  influences  whitecap  coverage. 

Monahan  et  al.  (1983)  suggest  that  whitecap  cover  depends  on  wind  dura¬ 
tion  and  fetch  as  well  as  on  speed  and,  in  shallow  water,  depends  on  the 
bathymetry  of  the  ocean  floor.  In  viscous  cold  water,  whitecap  coverage 
will  be  less  than  in  less  viscous  warm  water  under  the  same  meteorological 
conditions  because  the  ocean  is  more  efficient  in  dissipating  the  energy 
input  by  the  wind;  hence,  fewer  breaking  waves  are  created.  Wu  (1988) 
suggests  that  fwc  is  related  to  the  cube  of  the  friction  velocity  u*,  which  is 

proportional  to  the  rate  of  energy  supplied  by  the  wind  per  unit  area  of  the 
ocean  surface.  Using  the  friction  velocity  rather  than  the  wind  speed  im¬ 
plicitly  incorporates  atmospheric  stability  in  the  determination  of fwc- 
Piazzola  et  al.  (2002)  use  physical  relationships  to  incorporate  fetch  limi¬ 
tations  in  their  whitecap  cover  equation. 

Spray  concentration  profile 

de  Leeuw  et  al.  (in  preparation)  develop  a  new  theoretical  model  for  the 
relationship  between  dF/dr  and  the  droplet  concentration  distribution 
c(z,r)  as  a  function  of  height  z  and  droplet  radius  r.  This  model  builds  on 
earlier  work  by  Fairall  and  Larsen  (1984)  and  Hoppel  et  al.  (2002,  2005) 
and  uses  concepts  that  treat  particle  deposition,  as  originally  formulated 
by  Slinn  et  al.  (1978)  and  Slinn  and  Slinn  (1980).  We  sketch  just  the  rudi¬ 
ments  of  that  model  here. 
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Conceptually,  we  assume  that  all  spray  droplets  are  formed  at  some  refer¬ 
ence  height  h  above  mean  sea  level.  Because  spray  generation  is  closely  as¬ 
sociated  with  breaking  waves  and  the  mechanical  tearing  of  wave  crests,  a 
good  choice  for  h  is  0.5H1/3,  half  the  significant  wave  height.  We  assume 
that  for  o  <  z  <h 


dF_ 

dr 


—  Vdhch  ’ 


(10) 


where  Ft  is  the  net  turbulent  flux  of  spray  droplets.  Specifically, 


Ft  =  w'c'(r ) 


(ll) 


is  a  typical  vertical  Reynolds  flux,  where  w'  is  the  turbulent  fluctuation  in 
vertical  velocity  and  c'(r)  is  the  turbulent  fluctuation  in  the  concentration 
distribution  of  droplets  with  radius  r.  The  overbar  denotes  a  time  average 
that  typically  is  30-60  minutes.  As  with  most  Reynolds  fluxes  in  the  near¬ 
surface  atmosphere,  Ft  is  assumed  to  be  constant  with  height.  For  simplic¬ 
ity  c{h,r),  the  spray  concentration  distribution  at  height  h,  is  written  as  Ch. 
Also  in  Equation  10, 


vdh  =  - 


V  +v 

g  1  a 


1  +  : 


Ysl 

V. 


/J 


(12) 


is  the  deposition  velocity  for  droplets  of  radius  r  at  height  h.  Vg  is  the  ter¬ 
minal  fall  speed  of  droplets  with  radius  r,  always  assumed  positive,  Va  is  a 
diffusive  velocity  described  in  Appendix  A,  and 

fu  =  (h/iTr,,t"- ■  (13) 


Here  k  =  0.40  is  the  von  Karman  constant;  u.  is  the  friction  velocity,  a  dy¬ 
namic  velocity  scale  of  the  atmospheric  surface  layer;  and 


6  = 


Ay 


u. 


(14) 
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is  the  approximate  thickness  of  the  molecular  diffusion  layer  at  the  ocean 
surface,  where  v  is  the  kinematic  viscosity  of  air  and  A  is  a  constant  that 
we  take  as  5  (de  Leeuw  et  al.  in  preparation). 


For  droplets  above  the  source,  dF/dr  =  o,  and  the  governing  equation  is 


F  =  -V 

r  Vg 


Cz  fzhCh 


zh 


(15) 


where  cz  is  c(z,r)  and 


f,H=(z/hyr,lh 


(16) 


On  substituting  Equations  12  and  15  in  Equation  10,  we  obtain 


1+T f^-fzi) 

V  9 


1_  fzh)dF  /  dr 


1  +  7 
v  g 


V. 


(17) 


where 


/z6=(z/8)-Vfcl\  (IB) 

Notice  that  Equation  17  relates  the  spray  concentration  distribution  cz  to 
the  spray  generation  function  dF /dr  if  we  know  Ch,  the  concentration  at 
height  h. 

Spray  concentration  distribution  at  the  surface 

We  test  three  assumptions  for  estimating  the  surface  concentration  distri¬ 
bution  Ch  for  use  in  Equation  17  from  the  known  dF/dr.  Hoppel  et  al. 
(2005)  and  de  Leeuw  et  al.  (in  preparation)  suggest  that  the  upward  flux  of 
droplets  must  be  the  same  as  the  downward  flux,  which  means  that  the 
deposition  velocity  of  the  droplets  should  be  used  to  determine  the  con¬ 
centration  from  the  flux: 


dF 


Vdhch 


dr 


(19) 
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According  to  Equation  to,  this  assumption  implies  Ft  =  o. 

A  second  estimate  of  Ch  comes  from  using  the  droplet  ejection  velocity  Vej 
as  suggested  by  Monahan  (1968).  This  is  the  initial  upward  velocity  of  the 
jet  droplets  that  are  created  by  collapsing  bubbles.  Blanchard  (1963)  ob¬ 
tained  the  ejection  velocity  for  droplets  with  radii  between  1.5  and  335  pm 
from  laboratory  measurements.  This  ejection  velocity  takes  the  place  of 
Vdh  in  Equation  10: 


dF_ 

dr 


Vejch- 


(20) 


A  third  method  for  estimating  Ch  is  based  on  postulating  that  the  friction 
velocity  u,  is  the  appropriate  scale  to  obtain  Ch  from  dF/ dr.  Spray  droplets 

are  freed  from  the  sea  surface  by  bursting  bubbles  or  by  some  other  insta¬ 
bility  mechanism.  Once  free,  they  naturally  tend  to  fall  back  into  the  sea. 
Only  if  turbulence  in  the  air  suspends  them  do  they  have  any  appreciable 
residence  time.  In  models,  such  turbulent  suspension  and  dispersion  is 
parameterized  in  terms  of  the  Langevin  equation  (e.g.,  Edson  and  Fairall 
1994;  Pattison  and  Belcher  1999)  in  which  the  rate  that  droplets  move  up 
and  down  in  response  to  the  turbulence  is  related  to  aw,  the  standard  de¬ 
viation  of  the  wind’s  vertical  velocity  component.  Because  ow  1.25m ,  u* 

is  a  reasonable  scale  for  relating  Ch  and  dF/dr.  Thus, 


dF 

dr 


uxh. 


(21) 


In  all  three  cases,  cz  is  obtained  from  Ch  using  Equation  17  and  the  result 
that  Ft  =  o. 


Deposition,  ejection,  and  friction  velocities  are  compared  in  Figure  6.  The 
ejection  and  deposition  velocities  differ  by  more  than  3  orders  of  magni¬ 
tude  for  small  droplets  but  are  nearly  the  same  for  droplet  radii  of  a  few 
hundred  microns.  Friction  velocities  tend  to  be  between  the  ejection  and 
deposition  velocities,  varying  with  wind  speed  and  atmospheric  stability 
rather  than  with  droplet  size. 

The  effects  of  these  three  assumptions  for  determining  Ch  from  dF/ dr  are 
significant.  Note  that  for  a  given  flux,  higher  droplet  velocities  imply 
smaller  droplet  concentrations.  Calculated  spray  liquid  water  content  pro- 
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files  are  shown  in  Figure  7a  for  U10  =  20  m  s1.  Liquid  water  contents  based 
on  the  deposition  velocity  assumption  are  much  higher  than  those  ob¬ 
tained  from  using  either  the  ejection  velocity,  which  results  in  extremely 
small  liquid  water  contents,  or  friction  velocity.  Figure  7b  shows  the  cumu¬ 
lative  contribution  to  the  liquid  water  content  at  5  and  50  m  above  sea 
level  (asl)  over  the  range  of  droplet  sizes.  Because  of  the  extremely  small 
deposition  velocities  of  small  droplets,  they  contribute  significantly  to  the 
liquid  water  content  at  5  m,  and  even  more  so  at  50  m,  where  the  larger 
droplets  are  less  numerous.  On  the  other  hand,  large  droplets  have  a  sig¬ 
nificant  contribution  to  the  liquid  water  content  based  on  ejection  velocity 
because  of  their  relatively  low  initial  velocity. 


Figure  6.  Variation  of  deposition  and  jet  droplet  ejection  velocities  with  spray  droplet  radius; 

variation  of  friction  velocity  with  U10. 

We  can  compare  our  calculated  Ch  values  to  measured  values.  Concentra¬ 
tion  distribution  profiles  have  been  published  by  Monahan  (1968)  and  de 
Leeuw  (1986a, b,  1987).  Measured  Ch  values  for  wind  speeds  of  12.4  m  s1  in 
the  North  Sea  (de  Leeuw  1987),  13  m  s_1  (de  Leeuw  1986a)  in  the  North  At- 
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lantic,  11  m  s_1  near  Aruba,  and  16  m  s1  in  Buzzards  Bay  (Monahan  1968) 
are  compared  to  calculated  values  for  a  wind  speed  of  13  m  s1  at  25  m  asl 
with  a  2.5-m  significant  wave  height  in  Figure  8.  In  these  comparisons,  all 
droplet  characteristics  are  evaluated  in  terms  of  rgo,  the  droplet  radius  at 
equilibrium  for  a  relative  humidity  of  80%.  That  is,  in  Equation  19,  for  ex¬ 
ample,  dF/ dr  =  dF/ dr  so-  All  the  calculated  Ch  values  agree  well  for  droplets 
with  radius  greater  than  about  100  pm,  but  in  this  range  they  are  smaller 
than  Monahan’s  measured  values.  For  droplet  radii  less  than  20  pm,  the 
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Figure  7.  Results  of  concentration  calculation  at  the  ocean  surface  for  Uio  =  20  ms1; 
comparison  of  deposition,  ejection,  and  friction  velocity  assumptions:  (a)  liquid  water  content 
profile;  contribution  of  droplets  to  the  liquid  water  content,  over  the  range  of  droplet  radii,  at 

(b)  5  m  asl  and  (c)  50  m  asl. 
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ejection  velocity  Ch  is  too  small  by  up  to  an  order  of  magnitude  compared 
to  de  Leeuw’s  measured  values.  The  deposition  velocity  Ch  has  roughly  the 
same  shape  as  de  Leeuw’s  data  for  small  droplets  but  is  too  large  by  one  or 
two  orders  of  magnitude.  The  friction  velocity  Ch  has  the  same  magnitude 
as  the  de  Leeuw  values  for  small  droplets,  but  the  curve  is  very  flat  in  com¬ 
parison,  implying  that  droplets  with  radii  less  than  about  20  pm  are 
equally  numerous. 

We  further  test  these  three  concentration  assumptions  in  Section  4  by  us¬ 
ing  them  to  determine  the  spray  cloud  liquid  water  content  distribution 
during  icing  events  for  two  drill  rigs  in  Alaskan  waters.  We  can  then  com¬ 
pare  those  rates  to  the  observed  icing  severity. 
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Figure  8.  Measured  sea  spray  concentration  densities  at  the  ocean  surface  compared  to 
calculated  values  based  on  the  three  droplet  velocity  assumptions. 
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4  Icing  Cases 

Data  are  available  from  only  two  semi-submersible  drilling  rigs  in  sea 
spray  icing  events.  Severe  sea  spray  icing  affected  the  Ocean  Bounty  and 
the  nearby  Dan  Prince  during  the  winter  of  1979-1980  in  Cook  Inlet, 
Alaska.  Six  distinct  events  occurred  during  the  121-day  period  that  the 
Ocean  Bounty  was  on  station.  So  much  ice  accumulated  on  the  superstruc¬ 
ture  in  one  event  that  drilling  mud  had  to  be  offloaded  to  keep  the  vessel 
from  losing  buoyancy.  These  severe  events  instigated  icing  studies  by  Min¬ 
erals  Management  Service,  including  one  led  by  Cold  Regions  Research 
and  Engineering  Laboratory’s  David  Minsk  (now  retired)  on  the  semi- 
submersible  vessel  Sedco  708  on  the  North  Aleutian  Shelf  during  the  fall 
and  winter  of  1982-1983.  These  two  data  sets  are  presented  and  discussed 
here,  with  model  estimates  of  the  spray  icing  intensity. 

Sedco  708 

The  semi-submersible  drilling  rig  Sedco  708  (Atlantic  Venture)  was  de¬ 
ployed  by  ARCO  Alaska  for  the  Continental  Offshore  Stratigraphic  Test  at 
56.274°N,  i6i.976°W  in  53  m  of  water  (Figure  9).  The  Sedco  708  main 
deck  was  about  128  ft  above  the  bottom  of  the  pontoons.  The  normal  drill- 


Figure  9.  Location  of  Ocean  Bounty  and  Sedco  708. 
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ing  draft  of  85  ft  puts  the  main  deck  at  about  43  ft  (13  m)  above  the  ocean 
surface.  An  ice  accretion  observation  program  and  the  evaluation  of  candi¬ 
date  ice-phobic  coatings  were  included  as  part  of  the  drilling  program.  The 
information  provided  here  on  the  observation  program  is  from  Minsk 
(1984a,  1985).  Three  Rosemount  model  872DC  ice  detectors  were 
mounted  on  the  180-ft-tall  derrick,  with  a  fourth  at  about  20  ft  above  deck 
on  the  railing  on  the  roof  of  the  diving  bell  storage  area.  In  addition,  six 
arrays  of  8-in. -long,  -diameter  cylinders  were  mounted  just  below 

the  main  deck  in  accessible  locations,  with  reasonable  exposure  to  sea 
spray  icing.  Five  panels  with  candidate  ice-phobic  coatings  applied  to  l-ft- 
square  steel  plates  were  also  mounted  below  the  main  deck. 

Only  two  significant  storms  occurred  while  the  drilling  rig  was  on  station: 
one  in  early  December  1982  that  was  not  documented,  and  one  from 
January  3  through  January  8, 1983.  Hourly  weather  and  sea  conditions 
during  that  storm  are  shown  in  Figure  10.  Wind  speeds  reached  as  high  as 
20  m/s  with  gusts  over  25  m/s,  concurrent  with  air  temperatures  as  cold 
as  -io°C.  Measured  significant  wave  heights  (bottom  panel)  were  almost 
six  meters  with  good  agreement  between  the  measured  and  modeled  val¬ 
ues  from  Equation  8.  The  water  temperature  decreased  from  3-9°C  to 
3-3°C  during  this  period. 

No  icing  was  measured  by  the  Rosemount  ice  detectors  on  the  derrick  or 
the  railing.  Up  to  about  5  in.  of  ice  accumulated  on  the  diagonal  cylindrical 
trusses  in  the  middle  of  the  below-decks  structure,  but,  surprisingly,  there 
was  no  ice  accumulation  at  the  same  level  on  the  windward  side  of  the 
platform.  The  lower  third  of  the  trusses  was  kept  clear  by  waves;  above 
that  level  the  ice  thickness  decreased  with  height  as  the  bottom  of  the  deck 
was  approached.  Up  to  about  1  in.  of  ice  accumulated  on  the  cylinder  ar¬ 
rays.  That  is  only  a  rough  estimate  because  the  arrays  were  inaccessible 
during  the  storm.  A  total  of  about  30  tons  of  ice  accumulation  was  esti¬ 
mated,  based  on  the  observed  distribution  of  ice  thicknesses  and  the 
amount  of  ballast  that  had  to  be  pumped  to  maintain  trim. 

To  model  the  accretion  of  spray  ice  during  this  storm,  a  number  of  as¬ 
sumptions  had  to  be  made.  Assuming  the  anemometer  is  mounted  at  20  m 
above  the  sea  surface  results  in  the  calculated  10-m  wind  speed  shown  in 
the  middle  panel  of  Figure  10.  An  unstable  atmosphere  is  assumed,  with 
the  air  heated  from  below  by  the  warmer  water.  Thus  the  vertical  mixing  of 
the  spray  droplets  is  greater  than  it  would  be  in  stable  conditions,  leading 
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to  relatively  small  vertical  droplet  concentration  gradients.  The  vertical 
variation  in  the  liquid  water  content  is  calculated  based  on  the  droplet 
concentration  distribution  calculated  from  the  droplet  concentration  at  z  = 
0.5//1/3.  That  in  turn  is  determined  from  the  spray  droplet  flux  at  the 
ocean  surface  divided  by  either  the  droplet  deposition  velocity,  the  droplet 
ejection  velocity,  or  the  friction  velocity,  as  discussed  in  Section  4.  The  ic¬ 
ing  rate  for  each  hour  is  determined  for  cylinder  diameters  ranging  from 
0.5  to  30  cm,  incorporating  the  effect  of  collision  efficiency.  These  icing 
rates  are  based  on  the  assumption  that  the  ice  freezes  instantaneously  on 
the  windward  side  of  the  cylinder,  with  the  ice-covered  cylinder  retaining 
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Figure  10.  Measured  air  temperature,  measured  wind  speed  and  calculated  10-m  wind 
speed,  and  measured  and  calculated  significant  wave  height  for  Sedco  708. 
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the  same  projected  area  perpendicular  to  the  wind  direction  throughout 
the  storm.  The  accumulation  of  ice  at  approximately  the  level  where  ice 
was  observed  on  the  large  cylindrical  diagonal  trusses  is  shown  in  Figure 
11a  for  all  three  droplet  velocity  assumptions.  The  accumulation  of  ice  at 
the  level  of  the  cylinder  arrays  and  just  above  the  main  deck  level  is  shown 
in  Figure  lib.  These  results  can  be  compared  to  the  5  in.  and  1  in.  of  ice 
that  was  estimated  on  the  trusses  and  the  cylinders,  respectively.  Keep  in 
mind  that  the  25-mm  increments  on  the  vertical  axes  in  Figure  11  are 
about  1  in. 

Ocean  Bounty 

Spray  icing  of  the  semi-submersible  exploratory  drilling  rig  Ocean  Bounty 
operated  by  Phillips  Petroleum  Company  during  the  winter  of  1979-1980 
is  described  in  Williams  (1981),  Nauman  (1984),  Nauman  and  Tyagi 
(1985),  and  Minsk  (1984b).  The  Ocean  Bounty  was  on  station  near  Kam- 
ishak  Bay  in  Lower  Cook  Inlet,  20  km  from  shore  in  160  m  of  water.  Spray 
icing  also  occurred  on  the  nearby  drill  rig  Dan  Prince,  however  weather 
information  is  available  only  for  the  Ocean  Bounty  location.  In  this  region 
wind  speeds  can  vary  dramatically  over  short  distances  because  of  funnel- 
ing  and  channeling  of  winds  through  gaps  in  the  nearby  mountains  so  the 
conditions  at  the  location  of  the  Dan  Prince  might  have  been  significantly 
different.  The  Ocean  Bounty  is  107-m  long,  81-m  wide,  with  the  main  deck 
16  m  above  the  ocean  surface.  Anemometer  height  was  276  ft  (84  m)  asl. 

From  September  24, 1979,  to  April  26, 1980,  it  experienced  21  days  of  sea 
spray  icing  on  the  superstructure  in  six  icing  storms  lasting  between  1  and 
7  days  each.  During  these  events  meteorological  and  oceanographic  data 
as  well  as  ice  accumulation  rates  were  recorded  by  OCEANROUTES  mete¬ 
orologist-observers  every  2  hours  from  6  am  to  6  pm  and  at  10  pm  and  2 
am  for  a  total  of  nine  observations  each  day.  Ice  accumulations  were  re¬ 
corded  as  light,  moderate,  heavy,  very  heavy,  or  extreme,  corresponding  to 
roughly  2,  3,  4,  6,  or  10  in.  of  ice  accumulated  per  day,  respectively.  The 
water  temperature  ranged  from  4-4°C  to  5.6°C. 

Williams  (1981)  describes  typical  storms  in  the  Cook  Inlet  area.  Strong 
west-northwesterly  orographic  winds  spilling  through  Kamishak  Gap 
across  lower  Cook  Inlet  occur  with  a  trough  of  low  pressure  oriented  paral¬ 
lel  to  the  mountain  range  over  the  Inlet.  These  conditions  occur  with  rising 
pressures  following  an  intense  low,  with  cold  air  pooling  over  the  Iliamna 
Basin.  Strong  northeasterlies  are  associated  with  high  pressure  over  interi- 
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Figure  11.  Modeled  ice  accretion  on  Sedco  708,  (a)  on  a  30-cm-diameter  cylinder  at  5  and  10 
m  (like  the  diagonal  trusses),  (b)  on  a  2-cm-diameter  cylinder  (like  the  cylinder  arrays)  at  10 

and  15  m. 


ERDC/CRREL  TR-09-3 


25 


or  Alaska  and  a  low  southwest  of  the  Inlet.  Strong  winds  from  the  east- 
southeast  are  associated  with  low  pressure  systems  moving  northward  into 
the  Gulf  of  Alaska  or  Bering  Sea.  While  the  Ocean  Bounty  was  on  station, 
icing  was  primarily  associated  with  west-northwesterly  winds. 

Wave  height,  wind  speed,  air  temperature,  relative  humidity,  and  ice  ac¬ 
cumulation  estimates  are  provided  in  tables  grouped  by  wind  direction  in 
Nauman  (1984).  Concurrent  weather,  wave,  and  icing  conditions  can  be 
recreated  from  these  tables,  assuming  that  the  data  are  in  the  same  order 
in  all  five  tables.  However,  the  chronological  order  has  been  lost,  so  these 
observations  cannot  be  associated  with  particular  storms.  Air  temperature 
data  are  plotted  in  Figure  12  in  the  top  panel,  with  measured  wind  speed 
and  the  calculated  10-m  wind  in  the  middle  panel,  and  measured  and  cal¬ 
culated  significant  wave  height  in  the  bottom  panel.  Lines  are  drawn 
through  the  data  points  that  are  grouped  by  wind  direction  only  to  make 
the  plots  easier  to  read.  It  was  significantly  colder  and  windier  in  these  ic¬ 
ing  events  than  in  the  Sedco  708  event.  The  significant  discrepancy  be¬ 
tween  the  measured  and  calculated  significant  wave  heights  here  is  in  con¬ 
trast  to  the  Sedco  708  location.  The  wave  heights  are  constrained  by  the 
short  fetch,  particularly  for  winds  with  a  westerly  component. 

Modeled  icing  rates  on  2-cm-diameter  cylinders  at  10  and  40  m  asl  for  the 
three  sea  surface  concentration  distribution  estimates  are  shown  in  Figure 
13  along  with  the  rough  estimates  of  the  daily  rate  of  ice  accumulation. 
Note  that  the  icing  rate  at  40  m  is  not  significantly  smaller  than  the  icing 
rate  at  10  m,  near  the  level  of  the  main  deck  for  any  of  three  droplet  veloc¬ 
ity  assumptions.  The  vertical  scale  for  the  top  three  panels  in  this  figure  is 
equivalent  to  the  scales  for  the  bottom  panel:  1  in.  of  ice  accumulation  per 
day  is  essentially  1  mm  per  hour.  The  icing  rate  based  on  deposition  veloc¬ 
ity  is  too  high  by  a  factor  of  10.  On  the  other  hand,  the  rate  based  on  ejec¬ 
tion  velocity  is  near  zero  in  many  hours  when  significant  icing  was  ob¬ 
served.  The  friction-velocity-based  icing  rate  appears  to  best  match  the 
observations. 

Ocean  Bounty  time  series 

Late  in  this  study,  MMS  provided  a  file  of  the  scanned  Ocean  Bounty 
handwritten  daily  Oceanographic  and  Meteorological  Data  Summary 
sheets.  Some  of  the  pages  were  too  light  to  read,  and  on  some  generally 
readable  pages  some  entries  were  illegible.  Summary  sheets  for  a  few  days 
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were  missing.  Within  these  constraints,  we  created  an  electronic  file  of  the 
data  for  days  with  at  least  one  recorded  air  temperature  at  32°F  or  below. 


10°20°30o40o280°  290°  300°  310°  320°330° 


Figure  12.  Ocean  Bounty  air  temperature,  wind  speed  (measured  and  corrected  to  10  m  asl), 
and  measured  and  calculated  significant  wave  height.  Vertical  lines  group  the  values  by  wind 

direction. 

Most  weather  and  sea  parameters  were  recorded  every  2  hours  from  0600 
to  1800.  In  addition,  maximum  wave  height,  wind  direction,  wind  and 
gust  speed,  air  temperature,  and  barometric  pressure  were  recorded  at 
0200  and  2200.  Water  temperature  was  measured  at  1400  almost  every 
day.  When  sea  spray  resulted  in  ice  accreting  on  the  Ocean  Bounty,  ob¬ 
servers  also  visually  estimated  the  severity  of  icing,  designating  it  as  light, 
moderate,  heavy,  very  heavy,  or  extreme.  The  corresponding  numerical 
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values  of  2,  3,  4,  6,  and  10,  roughly  interpreted  as  inches  of  ice  per  day, 
were  recorded  on  the  data  sheets  between  0600  and  1800.  The  relevant 
oceanographic,  meteorological,  and  icing  data  from  the  handwritten  sum¬ 
mary  sheets  are  provided  in  Appendix  B.  The  21  days  with  spray  icing 
(Nauman  1984)  are  in  red. 
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Figure  13.  Icing  rate  on  2-cm-diameter  cylinders  at  10  (top  line  in  each  pair)  and  40  m  (lower 
line  in  each  pair)  on  the  Ocean  Bounty  for  the  three  droplet  velocity  assumptions,  compared 
to  the  estimated  icing  rate.  Values  are  grouped  by  wind  direction. 
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We  again  modeled  spray  icing  as  described  in  Section  3.  We  handled  miss¬ 
ing  data  (typically  at  0200  and  2200)  by  filling  in  with  data  from  the  near¬ 
est  hour.  Results  are  shown  in  Figure  14.  The  measured  wind  speed  cor¬ 
rected  to  10  m  asl  is  shown  in  the  top  panel  along  with  the  observed 
significant  wave  height  in  red.  The  second  panel  displays  the  air  tempera¬ 
ture  (black)  and  water  temperature  (red).  The  observed  icing  severity  is 
shown  in  the  third  panel,  with  the  icing  severity  observation  converted 
into  the  rough  equivalent  in  inches  per  day,  which  is  essentially  the  same 
as  millimeters  per  hour.  Modeled  icing  rates  using  the  ejection,  deposition, 
and  friction  velocity  to  calculate  the  liquid  water  content  distribution  from 
the  flux  are  shown  in  the  bottom  three  panels.  As  seen  in  the  previous  sec¬ 
tion,  using  the  deposition  velocity  results  in  excessively  high  icing  rates. 
Ejection  and  friction  velocity  results,  on  the  other  hand,  are  in  the  right 
range. 

Some  of  the  modeled  spray  icing  events  were  not  observed  on  the  Ocean 
Bounty.  The  most  significant  was  at  the  beginning  of  December  in  1979.  A 
possible  explanation  for  the  lack  of  accreted  ice  is  the  relatively  warm 
ocean  temperatures  (6°C)  and  air  temperatures  (-7°C)  in  the  days  from 
1979.94  to  1979.95.  This  result  suggests  that  incorporating  a  heat  balance 
calculation  in  the  spray  icing  model  would  improve  its  performance. 
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Figure  14.  Ocean  Bounty  time  series  with  modeled  spray  icing  rates  at  10,  25,  and  50  m  asl. 


ERDC/CRREL  TR-09-3 


30 


5  Discussion 

Droplet  concentration  at  the  surface 

A  crucial  required  input  to  an  ice  accretion  model  for  sea  spray  icing  is  the 
variation  of  the  spray  droplet  concentration  with  height  above  the  ocean 
surface.  While  we  can  specify  the  spray  droplet  flux  at  the  ocean  surface, 
determining  the  spray  concentration  proved  to  be  surprisingly  difficult. 
The  problem  boils  down  to  answering  the  question  “How  fast  are  the  drop¬ 
lets  moving?” 

Our  first  attempt— following  the  lead  of  others  in  the  sea  spray  field— was 
to  say  that  at  equilibrium  the  upward  flux  of  droplets  from  the  surface 
must  equal  the  downward  flux;  then,  the  spray  droplet  deposition  velocity 
is  the  correct  speed  to  use.  However,  this  results  in  droplet  concentration 
densities  at  the  ocean  surface  that  are  an  order  of  magnitude  too  high 
compared  to  measured  values  in  the  film  droplet  size  range  and  most  of 
the  jet  droplet  range.  And,  in  fact,  the  upward  flux  of  droplets  in  these  size 
ranges  may  be  significantly  different  than  the  downward  flux  as  the  drop¬ 
lets  evaporate  from  their  radius  at  formation  to  a  radius  in  equilibrium  at 
the  relative  humidity  of  the  air.  As  they  evaporate,  their  settling  velocity 
decreases,  and  they  are  carried  more  easily  aloft  by  the  turbulence  in  the 
air.  A  particular  droplet  may  remain  aloft,  evaporating  to  a  salt  particle,  or 
may  fall  back  into  the  ocean  miles  away  from  where  it  was  formed,  with  a 
smaller  radius. 

Our  second  attempt  was  to  use  the  measured  ejection  velocity  of  jet  drop¬ 
lets  as  they  are  formed  by  collapsing  bubbles.  This  incorrectly  treats  film 
and  spume  droplets  as  jet  droplets  and  assumes  that  the  droplets  retain 
that  high  (except  for  the  very  largest  droplets)  initial  velocity  over  a  sig¬ 
nificant  distance  rather  than  being  slowed  by  drag.  For  droplets  with  a  ra¬ 
dius  of  up  to  about  15  pm,  this  assumption  results  in  droplet  concentration 
densities  at  the  ocean  surface  that  are  an  order  of  magnitude  too  low 
compared  to  measured  values.  However,  ejection  and  deposition  velocities 
are  roughly  the  same  for  large  droplets,  so  spume  droplet  concentration 
densities  are  equally  well  represented  by  either  of  these  two  velocity  as¬ 
sumptions.  It  should  be  noted,  however,  that  the  calculated  concentration 
distribution  in  this  size  range  is  smaller  than  measured  values,  implying 
that  these  droplets  are  moving  more  slowly  than  we  expect. 
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Finally,  we  used  the  friction  velocity  to  parameterize  the  vertical  velocity  of 
the  droplets.  This  choice  is  intrinsically  different  from  the  other  two  ap¬ 
proaches  in  that  the  droplet  motion  is  associated  with  the  wind  and  at¬ 
mospheric  stability  rather  than  being  a  property  of  the  droplet.  This  choice 
provided  the  best  agreement  between  measured  and  calculated  concentra¬ 
tion  distributions  but,  like  the  other  two,  resulted  in  smaller  than  meas¬ 
ured  values  in  the  spume  range. 

The  friction  velocity  assumption  provided  the  best  agreement  with  the  lim¬ 
ited  information  we  had  for  two  cases  of  spray  icing  on  drill  rigs.  This 
choice,  however,  does  not  provide  the  observed  shape  of  the  ocean-level 
concentration  distribution.  Because  friction  velocity  does  not  vary  with 
droplet  size,  the  concentration  distribution  has  the  same  shape  as  the 
spray  generation  function  dF /  dr  &  0. 

Whitecaps  and  significant  wave  height 

Both  the  whitecap  coverage  and  the  significant  wave  height  can  be  calcu¬ 
lated  from  the  wind  speed  at  10  m.  The  parameterization  for  the  signifi¬ 
cant  wave  height  (Equation  8)  is  based  on  buoy  data  from  the  North  Atlan¬ 
tic.  This  formulation  provided  an  excellent  match  to  the  Sedco  708 
observations,  indicating  that  the  very  poor  match  to  the  Ocean  Bounty  ob¬ 
servations  is  likely  not  because  of  some  intrinsic  difference  between  wind 
waves  in  Alaskan  waters  and  in  the  Atlantic.  The  Ocean  Bounty  significant 
wave  height  discrepancy  probably  results  from  the  limited  fetch  at  the  drill 
rig’s  location  in  Lower  Cook  Inlet. 

The  relatively  small  significant  wave  heights  in  Cook  Inlet  might  also  im¬ 
ply  that  the  calculated  whitecap  coverage  from  Equation  9  is  higher  than 
actually  occurred.  The  maximum  Ocean  Bounty  wind  speeds  resulted  in 
calculated  whitecap  covers  of  nearly  100%.  The  relatively  small  wind 
waves  might  imply  fewer  breaking  waves  and  therefore  fewer  bubbles  and 
fewer  film  and  jet  droplets  than  is  expected  for  those  wind  speeds  in  open 
ocean  conditions.  However,  the  shearing  of  spume  droplets  off  wave  crests 
by  the  wind  should  not  be  diminished.  It  would  be  useful  to  test  the 
Piazzola  et  al.  (2002)  fetch-limited  whitecap  parameterization,  not  only 
with  the  Ocean  Bounty  data  but  also  with  whitecap  observations  in  limited 
fetch  conditions. 
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Wind  and  icing  data 

The  available  weather  and  icing  data  for  sea  spray  icing  on  offshore  struc¬ 
tures  was  very  limited.  We  still  expect  to  get  chronological  Ocean  Bounty 
weather  data  for  the  winter  of  1979-1980  from  MMS.  It  would  be  useful  to 
process  the  data  for  that  entire  winter  to  determine  how  well  our  algo¬ 
rithms  identify  the  significant  spray  icing  events  that  Nauman  (1984) 
listed. 

Available  wind  and  temperature  data  for  Alaskan  waters  could  be  used  to 
model  spray  ice  accretion  at  those  locations.  In  a  study  completed  recently 
for  MMS,  Veltcamp  and  Wilcox  (2007)  compiled  meteorological  data  from 
January  2001  through  September  2006  from  five  MMS  stations  along  the 
coast  of  the  Beaufort  Sea.  The  parameters  measured  included  wind  speed 
and  direction,  air  temperature,  humidity,  and  barometric  pressure.  Three 
of  the  stations  are  on  islands,  one  is  on  a  point  of  land  surrounded  by  wa¬ 
ter  on  three  sides,  and  one  is  2  km  inland  from  the  Beaufort  Sea.  All  except 
one  of  those  stations  is  still  operating  under  the  auspices  of  the  University 
of  Alaska.  Veltcamp  and  Wilcox  also  compiled  wind  data  from  29  addi¬ 
tional  stations  in  the  area,  some  of  which  were  in  operation  for  only  a  few 
months.  Some  of  these  hourly  data  extend  as  far  back  as  1984.  The  pur¬ 
pose  of  the  study  was  to  compile  weather  data  to  use  in  MMS’s  Oil  Spill 
Risk  Analysis,  Coastal  Zone  Oil  Spill,  and  oil  weathering  and  near-shore 
circulation  models.  Those  data  could  also  be  used  with  the  National  Snow 
and  Ice  Data  Center  maps  of  sea  ice  extent 

f http://nsidc.org/data/seaice  index/archives/index.htmD  and  this  sea  spray  icing  model  to 
determine  the  severity  of  spray  icing  on  offshore  structures  in  this  region. 
The  sea  ice  maps  would  be  used  to  screen  the  high  wind/cold  weather 
events  to  extract  only  those  that  occurred  with  open  water.  It  would  be  in¬ 
teresting  to  examine  the  history  of  modeled  spray  ice  thicknesses  to  iden¬ 
tify  any  changes  associated  with  global  warming. 
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Appendix  A 

Diffusive  velocity  Va 

From  de  Leeuw  et  al.  (in  preparation), 

r1/2  ici  i 

Va=^-io -3/stu.+  — . 

k  Rp 

Cdz  is  the  drag  coefficient  at  height  z.  St  is  the  Stokes  number 

Vu! 

St  =  -^— , 

vg 


(A-l) 


(A-2) 


where  g  is  the  acceleration  of  gravity.  Rp  is  a  molecular  diffusion  term 


Rp  =  In 


2ASc 


1/2 


IT 


+  ASc 


1/2 


(A-3) 


where  Sc  is  the  droplet  Schmidt  number 

„  v 

SC  =  —  .  (A-4) 


Dd  is  the  molecular  diffusion  of  droplets  (Pruppacher  and  Klett  1978,  p. 

361) 


D 


d 


kT(i  +  oJCn) 
6t\p  vr 


(A-5) 


where  k  ( =  1 .3806503  x  10  23  J/ K )  is  the  Boltzmann  constant,  T is  the  abso¬ 
lute  temperature  (in  Kelvin),  p  is  the  air  density,  and  r  is  the  droplet  ra¬ 
dius.  Kn  is  the  Knudsen  number 


Kn  = 


r 


(A-6) 
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where  Xa  is  the  mean  free  path  of  air  molecules.  Also  in  Equation  A-5 
(Pruppacher  and  Klett  1978,  p.  361), 


a  =  1.257  +  o.400exp 


-1.10 


Kn 


(A-7) 


We  calculate  these  parameters  as  well  as  u,  and  several  other  micrometeo- 
rological  quantities  in  the  marine  atmospheric  surface  layer  from  the  bulk 
flux  algorithm  developed  by  Andreas  et  al.  (2008)  using  the  mean  wind 
speed,  air  temperature  and  humidity,  surface  temperature  and  salinity, 
and  barometric  pressure. 
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Appendix  B 

Ocean  Bounty  time  series 


Date 

Hour 

SigHft 

Max  H  ft 

DirN 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

Icing 

11/19/79 

2 

20 

50 

58 

74 

34 

977.3 

6 

20 

32 

30 

75 

93 

33 

32 

968.5 

8 

22 

34 

40 

65 

90 

33 

32 

966.5 

10 

28 

41 

50 

58 

75 

33 

31 

964.4 

12 

23 

32 

60 

38 

50 

34 

33 

963.1 

14 

16 

24 

110 

32 

42 

40 

45 

36 

961.7 

16 

14 

22 

290 

25 

35 

36 

34 

962.1 

18 

11 

20 
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44 
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32 

31 
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22 

15 

300 

40 

47 

31 

965.5 
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7 

13 
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32 
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8 

6 

11 
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16 

24 

29 

27 
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10 

6 
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48 
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12 

5 
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24 
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31 

30 

983.4 

14 

4 

7 
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22 

30 

34 

44 

22 

984.8 

16 

4 

7 
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22 

30 

35 

30 

986.1 

18 

3 
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22 

28 

35 

24 

987.8 

22 

5 
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17 

22 

35 

989.2 

11/21/79 
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16 
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6 

11 

40 

47 
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38 

27 
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30 
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62 
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69 
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28 
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10 
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60 

70 

32 

45 

28 

966.5 

18 

10 

19 

290 

62 

69 

31 

30 

965.8 

22 

300 

58 

66 

31 

967.8 

11/22/79 

2 

300 

64 

72 

27 
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47 
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18 
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11 
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290 

46 

54 
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17 

978.7 

12 
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18 
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38 

46 

25 

19 

981.0 

14 

8 

15 
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37 
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25 

21 

983.1 

16 

7 

13 
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34 

44 

26 

45 

23 
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13 
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36 
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24 
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22 
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Date 

Hour 

SigHft 

Max  H  ft 

DirN 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

8 

4 

7 

270 

13 

19 

28 

26 

1000.7 

10 

4 

7 

240 

6 

29 

28 

1000.3 

12 

5 

8 

210 

19 

26 

31 

17 

1005.8 

14 

4 

7 

200 

9 

12 

31 

44 

20 

1006.8 

16 

5 

8 

280 

12 

17 

30 

23 

1008.1 

18 

4 

7 

290 

11 

17 

31 

26 

1009.8 

22 

300 

26 

32 

30 

1012.5 

11/24/79 

2 

310 

21 

27 

30 

1012.9 

6 

3 

5 

290 

19 

24 

29 

22 

1012.9 

8 

3 

4 

260 

6 

14 

30 

22 

1012.5 

10 

2 

3 

260 

5 

11 

31 

25 

1012.2 

12 

3 

4 

190 

17 

25 

31 

25 

1010.8 

14 

4 

6 

190 

29 

37 

32 

20 

1009.1 

16 

4 

7 

190 

34 

45 

32 

45 

16 

1008.1 

18 

5 

9 

170 

33 

42 

32 

19 

1006.1 

22 

160 

50 

59 

34 

1003.7 

11/25/79 

2 

140 

54 

61 

38 

999.7 

6 

9 

17 

130 

62 

72 

39 

28 

995.6 

8 

9 

17 

130 

61 

72 

39 

34 

995.3 

10 

8 

15 

130 

52 

62 

41 

35 

994.2 

12 

8 

16 

130 

54 

66 

42 

36 

992.9 

14 

8 

15 

130 

52 

64 

39 

45 

35 

991.2 

16 

9 

20 

280 

55 

62 

30 

29 

992.2 

18 

8 

16 

270 

54 

62 

31 

24 

993.2 

22 

52 

68 

30 

997.0 

12/07/79 

2 

30 

25 

40 

34 

1011.2 

6 

4 

7 

20 

24 

34 

33 

1009.8 

8 

4 

7 

340 

23 

34 

32 

1010.5 

10 

5 

8 

320 

27 

3 

31 

1010.2 

12 

5 

8 

310 

26 

38 

30 

1009.8 

14 

4 

7 

300 

25 

32 

31 

45 

20 

1008.8 

16 

3 

5 

310 

21 

28 

30 

22 

1008.5 

18 

4 

7 

300 

28 

36 

30 

22 

1008.5 

22 

4 

300 

8 

16 

30 

1007.8 

12/08/79 

2 

8 

290 

48 

54 

22 

6 

8 

15 

290 

53 

65 

21 

18 

1009.1 

8 

9 

17 

290 

55 

67 

21 

16 

1010.2 

10 

11 

20 

280 

54 

71 

21 

16 

1011.9 

12 

12 

23 

290 

69 

84 

21 

16 

1011.9 

14 

12 

23 

290 

71 

84 

22 

43 

16 

1014.2 

16 

13 

25 

290 

66 

81 

20 

16 

1016.9 

18 

14 

27 

290 

72 

83 

20 

16 

1018.6 

22 

25 

290 

58 

71 

1023.7 

12/09/79 

2 

13 

290 

46 

54 

24 

1024.4 

ERDC/CRREL  TR-09-3 


Date 

Hour 

SigHft 

Max  H  ft 

Dir  N 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

12/14/79 

2 

6 

360 

20 

29 

24 

1027.1 

6 

6 

11 

10 

29 

43 

24 

17 

1027.4 

8 

6 

11 

20 

39 

52 

25 

20 

1028.1 

10 

6 

11 

30 

34 

48 

25 

21 

1029.1 

12 

6 

11 

30 

34 

45 

29 

41 

22 

1030.1 

14 

6 

11 

30 

35 

45 

29 

25 

1030.8 

16 

6 

11 

30 

31 

44 

30 

24 

1031.2 

18 

6 

11 

40 

27 

46 

31 

28 

1033.2 

22 

10 

40 

24 

58 

1031.8 

12/15/79 

2 

10 

40 

26 

37 

32 

1031.2 

6 

6 

11 

40 

30 

42 

31 

28 

1029.5 

8 

6 

11 

40 

31 

41 

31 

28 

1028.4 

10 

6 

11 

40 

33 

44 

32 

30 

1027.1 

12 

6 

11 

40 

33 

45 

32 

28 

1025.4 

14 

7 

12 

40 

32 

48 

31 

42 

29 

1023.4 

16 

7 

13 

40 

42 

57 

30 

29 

1021.0 

18 

8 

15 

40 

43 

56 

29 

27 

1018.6 

22 

16 

20 

46 

58 

28 

26 

1013.2 

12/16/79 

2 

18 

20 

46 

58 

24 

1005.1 

6 

9 

17 

30 

45 

60 

27 

24 

999.7 

8 

10 

19 

40 

46 

67 

29 

24 

997.0 

10 

11 

21 

30 

52 

68 

30 

26 

993.9 

12 

11 

21 

20 

59 

75 

31 

27 

990.9 

14 

12 

22 

20 

62 

80 

31 

43 

25 

988.2 

16 

13 

25 

20 

65 

79 

31 

27 

986.1 

18 

13 

25 

20 

66 

80 

30 

27 

984.4 

22 

26 

20 

64 

81 

29 

981.0 

12/17/79 

2 

21 

20 

60 

79 

28 

978.3 

6 

13 

25 

10 

48 

69 

28 

25 

976.0 

8 

14 

27 

360 

42 

60 

28 

24 

975.3 

10 

12 

23 

360 

40 

56 

29 

25 

974.3 

12 

12 

23 

350 

42 

61 

26 

973.2 

14 

12 

23 

350 

40 

54 

28 

972.6 

16 

12 

23 

300 

54 

65 

20 

970.9 

18 

12 

23 

300 

62 

72 

17 

970.5 

22 

25 

300 

69 

78 

12 

970.2 

12/18/79 

2 

256 

300 

72 

82 

8 

970.9 

6 

15 

27 

300 

70 

83 

7 

972.6 

8 

15 

27 

290 

66 

80 

7 

973.9 

10 

15 

28 

290 

69 

82 

7 

975.3 

12 

15 

28 

290 

74 

93 

7 

974.9 

14 

16 

30 

290 

73 

92 

8 

977.3 

16 

16 

30 

290 

68 

83 

10 

978.3 

18 

16 

30 

290 

74 

86 

9 

978.7 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 


I-09-3 


42 


SigHft 

Max  H  ft 

Dir  N 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

Icing 

30 

290 

61 

73 

7 

979.0 

20 

270 

63 

72 

7 

979.3 

12 

23 

300 

53 

62 

9 

979.3 

12 

23 

300 

51 

59 

12 

979.6 

11 

20 

300 

44 

51 

13 

980.0 

8 

14 

300 

43 

51 

13 

980.4 

7 

12 

300 

43 

52 

12 

42 

980.4 

7 

12 

290 

44 

57 

12 

980.7 

7 

13 

290 

53 

67 

11 

980.7 

300 

65 

79 

8 

979.3 

300 

12 

23 

310 

78 

88 

8 

12 

24 

290 

79 

95 

9 

6 

14 

26 

290 

81 

100 

8 

6 

15 

28 

290 

88 

110 

7 

10 

15 

28 

290 

90 

109 

6 

10 

16 

30 

300 

95 

111 

5 

10 

16 

31 

290 

98 

114 

6 

10 

290 

92 

109 

10 

300 

81 

96 

13 

25 

290 

78 

85 

7 

6 

12 

22 

290 

76 

81 

6 

6 

10 

18 

290 

60 

77 

7 

6 

10 

18 

290 

63 

75 

6 

6 

10 

18 

290 

67 

78 

7 

42 

6 

10 

18 

290 

63 

75 

7 

6 

10 

17 

290 

61 

72 

8 

6 

300 

57 

64 

8 

300 

57 

70 

10 

17 

300 

60 

71 

5 

6 

9 

17 

290 

58 

67 

5 

6 

10 

18 

300 

64 

73 

6 

6 

10 

19 

290 

60 

70 

6 

6 

11 

19 

290 

64 

75 

5 

42 

10 

11 

19 

290 

64 

73 

3 

10 

11 

20 

290 

62 

75 

3 

10 

300 

57 

66 

3 

1004.9 

300 

51 

58 

3 

1005.4 

8 

15 

300 

43 

49 

7 

1004.1 

2 

8 

14 

290 

39 

47 

8 

1003.1 

2 

7 

13 

290 

32 

39 

13 

1003.1 

2 

6 

11 

300 

25 

32 

16 

1001.0 

5 

9 

360 

16 

25 

19 

42 

999.0 

4 

7 

360 

15 

23 

22 

998.6 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 


I-09-3 


43 


SigHft 

Max  H  ft 

Dir  N 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F  Tdew  °F 

P  mb 

Icing 

3 

5 

20 

15 

26 

22 

997.0 

40 

8 

21 

995.3 

140 

1 

22 

992.2 

3 

5 

40 

25 

28 

985.8 

5 

9 

50 

32 

44 

28 

982.4 

7 

13 

50 

51 

66 

27 

977.3 

12 

22 

120 

52 

57 

36 

971.9 

13 

25 

110 

64 

81 

36 

41 

966.8 

14 

27 

100 

74 

96 

35 

961.4 

17 

35 

100 

85 

113 

35 

958.0 

40 

51 

75 

34 

954.3 

30 

62 

85 

34 

950.9 

13 

24 

30 

50 

71 

29 

949.9 

10 

19 

350 

43 

58 

30 

951.9 

9 

17 

360 

27 

43 

31 

954.9 

8 

15 

50 

15 

22 

34 

957.3 

7 

13 

100 

17 

24 

33 

40 

958.0 

7 

13 

190 

28 

36 

34 

959.0 

6 

11 

180 

31 

36 

34 

961.4 

210 

14 

18 

33 

965.1 

160 

12 

17 

33 

964.5 

3 

5 

160 

8 

13 

36 

964.9 

3 

5 

320 

9 

15 

34 

970.9 

4 

7 

320 

18 

23 

33 

972.9 

5 

8 

300 

27 

32 

32 

973.6 

5 

9 

290 

31 

36 

32 

41 

975.3 

6 

10 

290 

34 

38 

30 

977.7 

5 

9 

290 

26 

30 

30 

979.0 

300 

24 

30 

29 

983.4 

300 

16 

23 

29 

986.1 

4 

6 

300 

20 

26 

26 

987.1 

4 

7 

310 

17 

22 

25 

988.5 

4 

7 

320 

20 

25 

25 

989.8 

4 

7 

310 

26 

32 

26 

990.2 

4 

7 

320 

28 

34 

26 

40 

990.9 

5 

9 

310 

35 

40 

26 

990.9 

5 

9 

310 

43 

49 

25 

990.5 

57 

64 

25 

990.2 

300 

57 

66 

23 

994.9 

9 

18 

290 

57 

67 

15 

997.6 

2 

12 

23 

290 

61 

72 

13 

999.7 

4 

13 

25 

290 

63 

75 

11 

1002.0 

6 

15 

28 

290 

64 

74 

10 

1004.1 

10 

14 

26 

290 

64 

74 

10 

40 

1006.1 

10 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 


I-09-3 


44 


SigHft 

Max  H  ft 

Dir  N 

Wind  kt 

Gust  kt 

Tair°F 

Twater°F  Tdew°F  P  mb 

Icing 

12 

23 

290 

61 

71 

10 

1009.5 

10 

12 

23 

280 

57 

72 

7 

1011.5 

10 

290 

57 

68 

6 

1016.9 

300 

55 

64 

5 

9 

17 

300 

55 

65 

6 

6 

9 

17 

290 

56 

64 

7 

6 

8 

15 

300 

50 

59 

7 

6 

7 

13 

300 

54 

62 

8 

6 

7 

13 

300 

50 

61 

8 

41 

6 

7 

13 

300 

52 

60 

7 

6 

7 

13 

290 

51 

59 

7 

1025.7 

6 

300 

48 

56 

6 

1025.4 

300 

44 

49 

8 

1025.1 

7 

13 

300 

46 

52 

10 

1023.7 

6 

7 

13 

300 

44 

50 

10 

1022.4 

6 

7 

12 

290 

45 

52 

9 

1021.7 

6 

7 

13 

290 

44 

51 

11 

1020.7 

6 

7 

13 

300 

43 

49 

10 

40  1019.3 

6 

7 

13 

300 

43 

49 

10 

1019.0 

6 

6 

11 

300 

40 

47 

11 

1019.0 

4 

300 

42 

51 

10 

1018.0 

300 

46 

65 

12 

1017.6 

7 

13 

310 

44 

52 

13 

1016.6 

7 

13 

310 

39 

47 

14 

1016.9 

7 

13 

320 

38 

47 

15 

1017.3 

7 

12 

320 

34 

46 

15 

1017.3 

6 

11 

360 

26 

31 

19 

41  1019.3 

6 

11 

360 

30 

30 

19 

1017.3 

5 

9 

360 

32 

38 

18 

1017.3 

360 

36 

54 

19 

1018.3 

360 

33 

48 

20 

1018.0 

6 

11 

10 

34 

44 

20 

1016.9 

6 

11 

360 

33 

44 

20 

1016.0 

5 

10 

360 

31 

43 

22 

1016.6 

6 

11 

350 

30 

40 

23 

1015.1 

6 

11 

360 

38 

46 

22 

41  1015.6 

6 

12 

350 

38 

49 

23 

1015.6 

6 

11 

320 

36 

46 

22 

1015.6 

330 

29 

37 

24 

1015.9 

310 

40 

46 

21 

1014.9 

5 

8 

310 

42 

48 

19 

1014.6 

5 

9 

310 

43 

50 

17 

1014.6 

6 

11 

300 

45 

55 

19 

1015.6 

6 

11 

310 

37 

44 

21 

1015.6 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 


I-09-3 


45 


SigHft 

Max  H  ft 

DirN 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

Icing 

6 

11 

310 

33 

38 

21 

41 

1015.2 

6 

11 

300 

39 

46 

18 

1016.6 

6 

11 

300 

40 

45 

16 

1017.3 

300 

40 

44 

20 

1019.0 

320 

32 

42 

21 

1021.7 

5 

9 

310 

31 

36 

21 

16 

1023.7 

5 

8 

320 

27 

33 

21 

17 

1024.7 

4 

6 

350 

20 

25 

25 

17 

1025.7 

4 

6 

360 

18 

26 

27 

13 

1026.8 

4 

6 

202 

21 

31 

27 

40 

18 

1026.8 

5 

8 

20 

26 

35 

26 

18 

1027.4 

5 

9 

30 

29 

42 

25 

16 

1027.8 

30 

37 

49 

24 

1028.4 

40 

34 

45 

30 

1028.1 

5 

9 

50 

30 

43 

31 

25 

1027.4 

6 

12 

130 

30 

31 

29 

1026.8 

6 

11 

110 

33 

38 

33 

31 

1026.8 

6 

12 

120 

39 

46 

35 

33 

1025.1 

7 

13 

120 

46 

55 

37 

40 

35 

1023.0 

8 

15 

110 

50 

58 

38 

36 

1021.3 

8 

15 

100 

48 

59 

39 

38 

1019.3 

110 

61 

72 

37 

1014.2 

320 

25 

29 

35 

1033.5 

3 

5 

300 

30 

36 

35 

22 

1028.1 

4 

7 

280 

32 

42 

35 

21 

1024.7 

6 

11 

290 

52 

64 

37 

23 

1021.0 

8 

15 

290 

63 

72 

36 

25 

1018.3 

8 

15 

290 

64 

77 

34 

39 

22 

1014.6 

9 

17 

300 

67 

77 

33 

18 

1012.5 

9 

17 

290 

66 

78 

27 

16 

1010.2 

320 

40 

57 

26 

1014.6 

20 

62 

84 

26 

1013.5 

9 

17 

10 

60 

90 

23 

11 

1013.9 

10 

19 

10 

58 

85 

23 

11 

1015.2 

4 

10 

19 

30 

61 

85 

22 

12 

1016.3 

4 

10 

18 

30 

56 

76 

21 

12 

1016.6 

4 

9 

17 

20 

47 

67 

20 

10 

1015.9 

4 

8 

15 

20 

48 

67 

19 

3 

1014.9 

4 

8 

15 

10 

44 

66 

19 

5 

1013.5 

4 

350 

40 

57 

20 

1011.2 

4 

310 

43 

52 

17 

1007.1 

8 

15 

310 

44 

53 

15 

10 

1002.4 

7 

14 

320 

39 

49 

14 

11 

1000.4 

6 

8 

15 

310 

48 

57 

12 

10 

998.0 

6 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 


I-09-3 


46 


SigHft 

Max  H  ft 

DirN 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F  Tdew  °F 

P  mb 

Icing 

8 

15 

310 

52 

64 

11 

9 

995.6 

6 

9 

17 

310 

56 

68 

13 

11 

992.2 

6 

9 

17 

300 

55 

68 

10 

8 

990.9 

6 

9 

17 

300 

55 

69 

9 

7 

987.5 

6 

340 

48 

78 

8 

983.4 

6 

300 

61 

74 

9 

979.0 

9 

17 

300 

62 

75 

6 

4 

976.0 

10 

18 

300 

66 

79 

6 

4 

974.6 

10 

10 

19 

290 

72 

89 

4 

2 

973.9 

10 

10 

18 

290 

69 

79 

4 

2 

973.2 

10 

10 

19 

290 

68 

81 

3 

1 

973.6 

10 

10 

19 

290 

64 

77 

3 

1 

975.3 

10 

10 

19 

300 

66 

70 

2 

-1 

976.0 

10 

10 

290 

61 

72 

2 

976.0 

10 

300 

58 

70 

1 

976.0 

9 

17 

300 

56 

67 

-1 

-3 

974.3 

9 

17 

300 

57 

66 

-2 

-4 

974.3 

10 

9 

17 

300 

58 

76 

-3 

-5 

972.9 

10 

9 

17 

290 

60 

74 

-2 

-4 

972.2 

10 

9 

17 

290 

62 

72 

-4 

-6 

971.2 

10 

10 

18 

290 

65 

76 

-5 

-7 

970.5 

10 

10 

18 

280 

63 

73 

-5 

-7 

969.9 

10 

290 

57 

68 

-3 

969.9 

10 

290 

54 

61 

-4 

970.9 

8 

15 

290 

51 

60 

-1 

-3 

969.9 

8 

15 

290 

51 

60 

-1 

-2 

969.9 

10 

8 

15 

290 

57 

66 

-1 

-2 

969.8 

10 

9 

17 

300 

65 

73 

0 

-1 

966.8 

10 

10 

19 

290 

75 

85 

2 

-1 

965.1 

10 

11 

21 

290 

72 

83 

1 

0 

966.8 

10 

11 

21 

300 

66 

78 

3 

1 

968.8 

10 

290 

64 

74 

4 

970.2 

10 

300 

68 

77 

972.2 

12 

22 

290 

80 

94 

-1 

-3 

974.9 

12 

22 

300 

79 

90 

-2 

-3 

977.7 

10 

12 

22 

290 

72 

81 

-4 

-6 

980.7 

10 

12 

23 

290 

78 

91 

-3 

-5 

981.0 

10 

13 

25 

290 

78 

90 

-4 

-6 

983.1 

10 

14 

27 

290 

78 

98 

-3 

-6 

984.4 

10 

13 

25 

300 

77 

86 

-4 

-5 

986.1 

10 

300 

65 

75 

-3 

990.9 

10 

300 

61 

68 

2 

992.9 

7 

13 

300 

35 

41 

7 

5 

994.9 

4 

6 

320 

9 

15 

10 

7 

994.9 

6 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 


I-09-3 


47 


SigHft 

Max  H  ft 

Dir  N 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

Icing 

4 

6 

170 

11 

17 

21 

18 

993.9 

3 

5 

270 

13 

21 

19 

993.6 

3 

5 

290 

4 

10 

18 

41 

17 

993.6 

4 

7 

310 

16 

24 

16 

15 

994.2 

5 

8 

320 

31 

36 

18 

17 

994.2 

320 

25 

30 

18 

995.6 

50 

20 

30 

30 

1037.6 

5 

9 

40 

30 

40 

29 

25 

1035.4 

6 

11 

40 

32 

41 

27 

26 

1037.3 

7 

14 

40 

29 

39 

28 

27 

1038.6 

7 

14 

40 

31 

41 

29 

27 

1038.6 

7 

13 

40 

32 

43 

28 

40 

27 

1038.6 

7 

12 

40 

33 

42 

29 

27 

1037.6 

7 

13 

40 

34 

44 

29 

27 

1037.3 

40 

34 

46 

29 

1036.9 

290 

47 

56 

26 

8 

15 

290 

47 

56 

23 

13 

994.3 

8 

15 

290 

52 

62 

23 

13 

997.3 

6 

9 

16 

290 

55 

66 

23 

12 

1001.7 

6 

11 

21 

290 

58 

68 

23 

14 

1004.4 

6 

11 

20 

290 

58 

67 

21 

14 

1008.1 

6 

9 

17 

300 

43 

50 

21 

40 

17 

1010.8 

6 

8 

15 

310 

39 

48 

19 

15 

1012.9 

4 

35 

42 

18 

1014.0 

4 

300 

29 

37 

19 

1035.1 

6 

11 

310 

33 

38 

12 

6 

10 

310 

28 

37 

12 

9 

2 

6 

10 

300 

28 

33 

20 

8 

1029.8 

2 

6 

11 

300 

24 

29 

21 

9 

2 

5 

9 

300 

21 

26 

22 

6 

1032.2 

2 

4 

6 

320 

7 

23 

40 

9 

1032.5 

3 

4 

350 

3 

7 

24 

17 

1033.5 

5 

11 

25 

1033.2 

50 

20 

30 

30 

25 

1037.6 

5 

9 

40 

30 

40 

29 

26 

1035.9 

6 

11 

40 

32 

41 

27 

27 

1037.3 

7 

14 

40 

29 

39 

28 

27 

1038.6 

7 

14 

40 

31 

41 

29 

27 

1038.6 

7 

13 

40 

32 

43 

28 

40 

27 

1038.6 

7 

12 

40 

33 

42 

29 

27 

1037.6 

7 

13 

40 

34 

44 

29 

27 

1037.3 

40 

39 

46 

29 

1036.9 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 


I-09-3 


48 


SigHft 

Max  H  ft 

DirN 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb  Icing 

130 

31 

37 

1036.9 

7 

13 

120 

22 

27 

36 

33 

1036.6 

6 

12 

110 

28 

33 

36 

35 

1036.6 

6 

11 

110 

24 

30 

37 

36 

1037.6 

5 

9 

110 

24 

30 

38 

35 

1037.9 

6 

10 

100 

25 

28 

38 

39 

35 

1037.3 

5 

9 

140 

25 

30 

39 

33 

1037.6 

5 

10 

140 

16 

39 

33 

1037.6 

140 

10 

20 

38 

1037.9 

320 

44 

54 

1026.8 

6 

11 

320 

56 

65 

31 

24 

1025.4 

6 

12 

320 

51 

61 

32 

25 

1025.1 

6 

11 

320 

56 

66 

32 

29 

1025.1 

6 

11 

320 

47 

56 

34 

26 

1024.0 

6 

10 

320 

47 

55 

35 

40 

26 

1023.0 

6 

10 

320 

48 

55 

33 

24 

1022.0 

6 

11 

310 

53 

59 

31 

25 

1021.3 

300 

56 

70 

33 

1015.9 

30 

53 

68 

32 

1015.2 

6 

11 

20 

48 

56 

32 

24 

1012.9 

7 

13 

30 

57 

67 

31 

22 

1010.8 

9 

17 

30 

69 

87 

31 

22 

1009.1 

10 

18 

40 

66 

82 

32 

26 

1009.1 

13 

25 

30 

82 

97 

31 

40 

26 

1007.1 

13 

26 

40 

75 

91 

28 

22 

1006.4 

14 

23 

40 

81 

94 

27 

20 

1005.8 

2 

40 

83 

94 

25 

1006.4 

2 

40 

80 

93 

25 

14 

28 

40 

85 

101 

26 

23 

1002.0 

4 

14 

27 

40 

85 

96 

25 

15 

1001.4 

4 

14 

27 

30 

82 

96 

26 

20 

1000.0 

4 

13 

25 

40 

77 

91 

26 

16 

999.7 

4 

13 

26 

40 

70 

81 

26 

40 

17 

998.0 

4 

13 

25 

40 

69 

85 

26 

20 

996.6 

4 

13 

26 

40 

73 

86 

23 

18 

995.9 

4 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 


I-09-3 


49 


SigHft 

Max  H  ft 

DirN 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

Icing 

30 

74 

91 

23 

989.5 

11 

21 

20 

67 

77 

24 

14 

986.5 

11 

20 

20 

64 

70 

25 

16 

985.1 

4 

10 

19 

330 

51 

63 

23 

17 

984.1 

4 

8 

15 

330 

50 

59 

21 

17 

982.4 

4 

8 

15 

330 

52 

70 

22 

41 

16 

980.7 

3 

9 

17 

320 

56 

69 

21 

17 

979.0 

3 

9 

17 

320 

60 

68 

21 

17 

978.0 

3 

320 

59 

69 

21 

976.0 

4 

310 

68 

76 

21 

973.9 

10 

19 

310 

63 

71 

18 

17 

973.9 

9 

17 

310 

59 

69 

18 

15 

974.9 

4 

9 

17 

310 

59 

66 

18 

17 

974.9 

4 

8 

15 

310 

49 

56 

18 

17 

974.9 

4 

8 

15 

320 

43 

54 

20 

41 

18 

973.2 

4 

8 

15 

310 

48 

53 

20 

18 

973.6 

4 

8 

15 

310 

41 

50 

20 

19 

974.9 

4 

320 

35 

41 

20 

976.0 

4 

300 

19 

29 

38 

971.2 

8 

15 

130 

48 

58 

38 

33 

968.5 

7 

13 

120 

38 

46 

39 

32 

968.2 

6 

11 

110 

32 

41 

39 

32 

967.5 

5 

9 

300 

21 

33 

30 

966.5 

5 

9 

40 

20 

31 

36 

41 

34 

964.8 

5 

9 

20 

22 

31 

34 

33 

964.4 

7 

13 

310 

50 

55 

34 

33 

964.4 

300 

51 

57 

33 

966.1 

300 

29 

36 

23 

970.5 

6 

11 

170 

37 

47 

27 

24 

976.3 

3 

5 

240 

15 

31 

30 

25 

977.3 

2 

4 

280 

8 

16 

28 

25 

978.3 

3 

5 

230 

16 

24 

28 

22 

979.7 

4 

7 

220 

31 

39 

28 

41 

21 

983.1 

3 

5 

220 

23 

31 

28 

19 

984.1 

3 

5 

190 

16 

26 

27 

23 

986.5 

230 

27 

44 

27 

992.2 

ERDC/CRREL  TR-09-3 


Date 

Hour 

SigHft 

Max  H  ft 

Dir  N 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

18 

22 

02/05/80 

2 

230 

45 

50 

37 

988.8 

6 

6 

11 

150 

37 

47 

32 

15 

991.9 

8 

3 

6 

310 

18 

23 

29 

20 

992.9 

10 

3 

5 

320 

19 

25 

29 

21 

992.9 

12 

4 

7 

320 

32 

37 

30 

24 

991.9 

14 

4 

7 

340 

27 

35 

30 

42 

27 

991.9 

16 

4 

7 

330 

35 

42 

28 

27 

988.8 

18 

5 

9 

340 

33 

41 

31 

30 

987.8 

22 

50 

41 

45 

29 

984.4 

02/06/80 

2 

50 

34 

38 

38 

980.0 

6 

8 

15 

120 

55 

63 

41 

40 

976.3 

8 

8 

15 

110 

50 

57 

41 

40 

976.6 

10 

8 

15 

130 

49 

58 

40 

39 

976.6 

12 

8 

15 

140 

50 

63 

40 

38 

977.0 

14 

9 

17 

140 

55 

64 

39 

43 

38 

977.3 

16 

11 

21 

230 

77 

88 

33 

29 

980.0 

18 

12 

22 

220 

79 

90 

30 

27 

983.7 

22 

280 

50 

65 

26 

02/07/80 

2 

300 

12 

26 

34 

998.3 

6 

7 

13 

230 

39 

49 

38 

34 

1004.4 

8 

3 

9 

220 

34 

41 

39 

31 

1006.4 

10 

4 

7 

180 

28 

36 

38 

26 

1009.1 

12 

4 

7 

190 

25 

37 

32 

26 

1010.5 

14 

4 

6 

220 

26 

30 

34 

42 

28 

1011.9 

16 

3 

5 

190 

20 

25 

34 

27 

1012.2 

18 

4 

7 

170 

31 

39 

34 

26 

1012.5 

22 

160 

45 

32 

36 

1013.2 

02/16/80 

2 

320 

14 

17 

34 

1016.3 

6 

2 

3 

310 

14 

17 

33 

28 

1016.6 

8 

1 

2 

330 

15 

18 

32 

30 

1017.6 

10 

2 

4 

340 

15 

18 

33 

28 

1013.9 

12 

2 

4 

350 

10 

13 

32 

26 

1012.5 

14 

2 

4 

330 

5 

33 

41 

25 

1011.2 

16 

2 

4 

310 

15 

18 

34 

26 

1010.2 

18 

2 

3 

310 

14 

19 

30 

26 

1009.5 

22 

330 

16 

25 

32 

1009.1 

02/17/80 

2 

6 

3 

4 

330 

29 

34 

28 

20 

8 

3 

5 

320 

33 

40 

28 

24 

10 

3 

4 

330 

21 

27 

28 

25 

12 

3 

5 

340 

13 

20 

29 

25 

999.3 

ERDC/CRREL  TR-09-3 


Date 

Hour 

SigHft 

Max  H  ft 

Dir  N 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

12 

5 

9 

50 

36 

46 

34 

29 

993.2 

14 

7 

15 

90 

36 

44 

36 

41 

33 

991.9 

16 

8 

15 

40 

54 

68 

34 

31 

988.5 

18 

9 

18 

60 

62 

73 

34 

32 

985.4 

22 

110 

78 

108 

34 

978.7 

02/23/80 

2 

130 

54 

64 

38 

979.0 

6 

13 

25 

130 

73 

90 

42 

38 

975.3 

8 

12 

23 

140 

72 

89 

41 

39 

977.3 

10 

11 

22 

140 

55 

68 

41 

37 

979.0 

12 

9 

18 

140 

58 

68 

41 

37 

980.0 

14 

9 

17 

150 

51 

60 

41 

40 

36 

983.1 

16 

9 

18 

220 

56 

69 

35 

27 

985.8 

18 

9 

19 

230 

67 

78 

34 

26 

991.9 

22 

10 

320 

12 

18 

30 

1001.0 

02/24/80 

2 

290 

17 

23 

29 

1006.8 

6 

3 

4 

350 

11 

16 

30 

26 

1008.1 

8 

4 

7 

10 

18 

23 

31 

26 

1008.1 

10 

5 

8 

40 

29 

35 

32 

31 

1008.1 

12 

5 

9 

40 

32 

38 

37 

35 

1007.5 

14 

7 

13 

110 

33 

38 

38 

40 

31 

1007.1 

16 

7 

13 

140 

20 

27 

37 

32 

1007.1 

18 

7 

13 

50 

39 

43 

35 

29 

1007.5 

22 

40 

31 

36 

35 

1008.1 

02/25/80 

2 

50 

24 

29 

36 

1010.5 

6 

4 

6 

40 

27 

30 

35 

28 

1010.2 

8 

5 

8 

40 

31 

36 

35 

30 

1009.2 

10 

5 

9 

30 

34 

39 

35 

28 

10096.1 

12 

6 

11 

30 

36 

42 

34 

28 

10096.1 

14 

6 

11 

20 

38 

43 

34 

41 

29 

1008.8 

16 

7 

13 

10 

41 

48 

34 

23 

1007.8 

18 

8 

15 

360 

49 

57 

34 

22 

1007.8 

22 

330 

31 

37 

36 

1006.1 

03/09/80 

2 

290 

48 

59 

34 

1007.8 

6 

7 

13 

290 

45 

50 

32 

23 

1008.1 

8 

7 

13 

290 

41 

49 

32 

25 

1008.8 

10 

7 

13 

300 

44 

52 

32 

25 

1009.1 

12 

7 

13 

290 

47 

53 

32 

25 

1009.1 

14 

7 

13 

300 

45 

52 

33 

41 

25 

1008.8 

16 

6 

11 

300 

40 

44 

33 

25 

1009.1 

18 

6 

11 

300 

44 

49 

32 

26 

1009.1 

22 

300 

42 

48 

32 

1009.1 

03/10/80 

2 

300 

36 

44 

34 

1009.5 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 

6 

8 

10 

12 

14 

16 

18 

22 

2 


!-09-3 


53 


SigHft 

Max  H  ft 

Dir  N 

Wind  kt 

Gust  kt 

Tair°F 

Twater  °F 

Tdew  °F 

P  mb 

5 

9 

300 

42 

47 

31 

23 

1007.8 

5 

9 

300 

43 

49 

32 

23 

1008.1 

5 

9 

310 

43 

47 

32 

24 

1008.1 

6 

11 

310 

42 

48 

32 

41 

24 

1007.8 

5 

9 

310 

40 

45 

33 

21 

1007.1 

5 

9 

300 

36 

43 

33 

24 

1006.4 

310 

38 

46 

31 

1005.8 

310 

41 

48 

30 

1005.1 

5 

9 

310 

38 

45 

30 

16 

1005.1 

5 

9 

310 

41 

43 

30 

17 

1005.1 

5 

9 

310 

40 

45 

31 

23 

1005.4 

6 

11 

310 

44 

46 

31 

23 

1005.8 

7 

13 

300 

49 

52 

30 

41 

25 

1005.8 

7 

13 

310 

54 

55 

29 

26 

1005.8 

8 

15 

300 

54 

62 

28 

22 

1005.8 

300 

55 

63 

26 

1006.4 

300 

58 

66 

25 

1007.8 

10 

19 

300 

67 

75 

23 

18 

1009.5 

11 

21 

300 

70 

79 

22 

17 

1009.8 

11 

21 

300 

69 

78 

23 

19 

1010.5 

11 

21 

300 

73 

72 

23 

17 

12 

22 

300 

73 

72 

23 

41 

18 

13 

25 

300 

71 

83 

22 

13 

1010.2 

13 

26 

300 

75 

85 

21 

13 

1010.2 

300 

74 

86 

19 

1010.9 

290 

76 

86 

15 

1009.5 

14 

27 

300 

85 

88 

15 

11 

1008.5 

14 

28 

300 

72 

91 

13 

10 

1008.8 

14 

28 

300 

77 

87 

15 

13 

1008.8 

13 

27 

290 

73 

89 

15 

40 

11 

10058.5 

13 

25 

290 

78 

91 

14 

10 

1008.5 

13 

25 

300 

86 

90 

14 

11 

1008.5 

13 

25 

300 

81 

91 

14 

10 

1008.5 

290 

70 

85 

13 

1008.5 

290 

66 

77 

13 

1007.5 

12 

22 

300 

62 

75 

12 

10 

1006.8 

7 

13 

300 
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